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INTRODUCTION 

The aim of this discussion is to reproduce the basic 
features of stellar structure and evolution (as found from 
accurate calculations) by purely analytic considerations 
in order to gain physical insight into the evolution of 
stars. We will not here attempt accurate calculations 
of structures and evolutionary tracks. First we discuss 
general properties of stellar structure and evolution. 

Then, analytic models. are constructed for the early homo- 
geneous and the advanced inhomogeneous stages of evolution. 
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I. EQUATIONS OF STELLAR STRUCTURE 

The basic equations governing the structure of stars 
are conservation of mass, conservation of momentum and 
conservation of energy ( Schwa rzschild, 1958 and Wrubel, 
1958). Rotation and magnetic fields will be neglected 
so that a star will be spherically symmetric. 

Hydrostatic Enuilibrium 


A star changes very slowly during most of its life 
and so may be considered in hydrostatic equilibrium. 

Two forces balance to keep a nonrotating star in hydrostatic 
equilibrium: the gravitational force directed inwards 

and the gas and radiation pressure force directed outward. 
The equation of hydrostatic equilibrium is 


dP _ GM r o 

dr " r 2 


( 1 . 1 ) 


The total pressure is the sum of gas and radiation pressure, 

p = p gas + p rad . | 

i 

Fpr an ideal gas 

p gas ~ ]Jjj” P t » ( 1 . 2 ) 

where H = 1.6 7 x 10“ ^g is the mass of a proton and |l is the 
mean molecular weight. 


P rad = 3 


(1.3) 
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M r is the mass interior to r; the equation of mass con- 
servation is 

M r “ J5 4iTr 2 pdr, (l.4a) 

or, in differential form, 

! 

! 

dM_ 0 1 

~ 4rrr*p ... (lV4b) 

Energy Conservation 

i 

The total energy of an element of material is 


E = u + n-+ K, (1.5) 

where U is the internal energy of the gas. Cl is gravita- 
tional potential energy, and K is the kinetic energy of 
large scale mass motion, which we are neglecting here. 
The internal energy of a gas plus radiation -is 

U = — ~ NkT + - aT 1 *, 

Y-l p 9 . 


where y is the ratio of specific heats (y = 5/3 for a mon- 
atomic ideal gas) • The sources and sinks of energy are 
(l) energy release due to nuclear reactions^ and ( 2 ) energy 
transport into and out of the element of material. 

Let £ be the net release of energy per gram per second, 
and F be the energy flux. The equation of conservation of 
energy is then 


dE 

dt 


dU dQ 
dt dt 


— divF 
P 


(per gram per second) . The change of gravitational potential 
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energy as 


dH — — dW — P d v = “'pZ' dP • 


Define the luminosity L r as the total net energy flux 
through a spherical shell of radius r, so that 


L r = 4rrr 2 F 


Then the equation of energy conservation is 


dLp -g 

dr 


4 " r X e + H - f ] . f 1 - 6 ) 


Energy Transport 


Energy is transported by radiation and convection, 

and by conduction when the electrons are degenerate. 


4 rr r2 


= F , + F 

rad r conv 


(1.7) 


In the interior of a star, where the radiation is 
almost isotropic, the momentum balance for radiation is 


dr c 4TTr 2 * 

where Pr = -j aT4 is the radiation pressure, (up) -1 is 
the photon mean free path, and c is the velocity of light. 
That is, the force due to, the gradient of the radiation 
pressure is equal to the momentum absorbed from the radiation 
beam in passing through matter. Thus, in the interior 
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of the star, the radiative energy flux is 

v 

3 


4acT 


dT 

dr * 


( 1 . 8 ) 


rad “ 3xp 

and the gradient necessary to drive the radiation flux is 

3 *o L, 


dt 

dr 


4ac 4rrr^ 


(1.9) 


The convective flux is (Spiegel, 1965) > crudely, the 
energy fluctuation (excess or deficiency) of an element of 
gas, times its velocity, averaged over a spherical surface 
in the star, 

F 


conv 


- P Cpw6 


( 1 . 10 ) 


where w is the radial velocity fluctuation and 0 the tem- 
perature fluctuations in the matter. The velocity and 
temperature excess or deficiency of a convective element 
depend on the superadiabatic gradient 

£ " ' 

Because convection is an extremely efficient energy trans- 
port mechanism, the superadiabatic gradient is very small 
and the temperature gradient will be very nearly equal 
to‘ the adiabatic gradient. 


fdT 

- \ 

i dr 

V- 

\ d rf ad. 


dT 

dr 


'dr\ = T_dP , (1.12) 

s drj a d r p dr 


where r is the effective ratio of specific heats, including 
ionization, dissociation and radiation. Near the surface. 
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where the photon mean free path is long, there is a leakage 
of heat by radiation from the convective elements and the 
convective temperature gradient is greater than the adiabatic 
gradient. 

Stellar Structure 

.* 

Order of magnitude estimates of the density, pressure, 
and temperature of a star can easily be made from the con- 
dition of hydrostatic equilibrium. The mean density is - - 

JT - 4% • U.U) 

In the equation of hydrostatic equilibrium ( 1 . 1 ), setting 

dP/dr ** (P c - p q )/ r » 

where P 0 is the central and P« the surface pressure, gives 

i p c m (1.14) 

c R tc 

since P Q « P c . Let 0 ** p gas/ p > the ratio of S as pressure 

f 

to total pressure and assume" that the material of the star 
is a perfect gas. Then the central temperature is obtained 
from the perfect gas law (equation 1.2), 

m " U0H Pc U8H Gjj 

• e * k £ . w k R 

The mean energy generation rate is 


( 1 . 15 ) 
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For the sun 

L *=* 3*89 x 10 33 ergs/sec, 

M » 1.99 x 10 33 g, (1.16) 

R “ 6.95 x 10 10 cm . 

-Thus the internal conditions of stars are of the order of 
magnitude • 

5 “ 1,41 (%) (f 2 ) g/c " 3 

P c « 1.1 x 10 16 dynes/cm 2 

2 > 3 llo 7 “|)(r) * K 
5 = 1-9 Or) eres/s ' sec 

as functions of the stars' mass, radius, and luminosity 
given in solar units.; 

For a more detailed account of the restrictions 
imposed by hydrostatic equilibrium on stellar structure 
see Chandrasekhar (1939) • 

This section is concluded by presenting a derivation 

of the expression for. the gravitational potential energy 

* 

of a sphere of uniform density. The gravitational potential 
energy is 

R 

n - £ Jo * d M(r) “ -3 J p dV, ( 1 . 18 ) 

where $ is the gravitational- potential. 
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For a sphere of uniform density the equation of hydro- 
static equilibrium (1,1) is 


1 dP „ d_ /P\ _ d$ 

{f dr dr [ py dr * 


Upon integrating, using the boundary condition that P/p -* 0 
at the surface, we get 

- $ + $ s = — and $ *» — ; 

P B R 


then 




“ d M(r) 


1 GM 2 

2 R 



Thus the gravitational potential energy of a sphere of 
uniform density is s* 


3 GM 2 
5 R 


(1.19) 


•The absolute value of the gravitational potential energy 
in an actual star will be somewhat larger, but of the same 
order of magnitude. 
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II. 

A star is a self-gravitating mass of gas in space. 

•The evolutionary trend of internal stellar conditions is 
determined by hydrostatic equilibrium and its radiation of . 
energy away into space. The life history of a star is 
the progressive concentration of its mass towards its 
center, pulled by its own gravitational field. This con- 
traction ^releases gravitational energy, heats up the gas, 
and, as the gas becomes hotter, thermonuclear reactions 
among various nuclei become possible. At certain temperatures 
the thermonuclear reactions can supply the energy losses, 

the gas and radiation pressure can support the star, 'and 

/*" 

the gravitational contraction is temporarily halted. < 

« ; 

A necessary condition for hydrostatic equilibrium 

i 

is the virial theorem for a self-gravitating mass (Chandrasekhar, 

1939), 

2 K+ft= 3 (y - 1) U + ft - 0 # (2.1) 

i 

Here K is the total thermal energy of the mass, U is itis 

I 

internal energy, . and ft is its gravitational potential energy. 

The virial theorem requires that the thermal energy of a 

i 

star equal half the absolute value of its gravitational 
potential energy (since ft is intrinsically negative). 

As a star contracts and releases gravitational energy, Q 
becomes more negative, and the thermal energy must 

increase. Half of the gravitational energy that is released . 
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is stored as thermal energy, increasing the temperature 
in the interior of the star, and half is radiated away* 
The mean relation of temperature to density can be. 
derived from the virial theorem. For a sphere of gas 
whose internal pressure is given by the perfect gas law 
with ratio of specific heats Y *® 5/3 » the virial theorem 
(2.1) becomes 


2 U + XI = 0 


( 2 . 2 ) 


For a uniform density distribution the gravitational energy 


is 


n ■- 


3 GM 2 

5 R 


.and the internal energy is 


»- 4 


(2.3) 


where M/jlH is the number of particles* Thus 


- ' 1 |iH M 

T = 5 k G R 


(2.4) 


The mean density (equation 1.13) is 


3 M 

TT = — z , 

- - 4 " R 3 * 

SO 

1 [ 4ttV 3 GuH 2/3 1/3 

T “ 1 It; ~ m ' p ' • (2-5) 

Thus the relation between temperature and density for 
stars with negligible radiation pressure is 

T - 4.1 x 10 6 ugp 3 pV3 . K . j(2.6) 
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Here T and p are the local temperature and density at any 
point in the star. 

The above temperature-density relation does not hold 
for those stars whose internal pressures are predominantly 
governed by the radiation pressure* Define the boundary 
line between stars whose internal conditions obey the 
perfect gas law and those whose internal conditions are 
regulated by radiation by an. equality of pressures for the 
two cases, i.e., \ 




i 



or 

T '= 2.55 x 10 7 p 1 / 3 . | 

* I 


The boundary corresponding to this condition occurs at 
5*5 For heavier stars radiation pressure, is predominant 

In such cases Y = 4/3’. and the virial theorem gives U m — SI* 
Thus , 

U - V a T 4 = - — ( 2 . 7 ) 

Expressing R in terms of the mean density (1.9) we obtain 
the temperature-density relation j 

T = 1.92 x 10 7 p 1 ^ 3 . ( 2 . 8 ) 

The temperature depends on density as before (to the l /3 
power) but the effect of mass is less pronounced. 

The temperature-density relations ( 2 . 6 ) and ( 2 . 8 ) 
describe. the dependence of the temperature on the density 


* 
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inside a star. The evolution of stars consists of pro- 
gressive gravitational contraction, increasing the central 
density and temperature according to 

T * 0 



interrupted at times by central nuclear burning. Some 
simplified evolutionary tracks for internal stellar con- 
ditions are shown in Figure 1. 


When the central density of a star gets very large, 
the matter may become degenerate and the equation of state 
thus changes. The boundary of degeneracy in terms of density 
and temperature has the asymptotic forms for low and high 
density (nonrelativistic and relativistic energies). 



low density. 


(2.9) 


high density. 


The full boundary curve is derived by Chandrasekhar (1939) « 
This boundary is also plotted in Figure' 1. 

Stars of mass less than about 1.3 % enter the de- 
generate region. For these stars the pressure due to 
degenerate electrons is so high that further compression 
is no longer possible. This is essentially the end point 
in the evolution of a star of small mass. The star be- • 
comes a white dwarf, achieving in this process some max- 
imum temperature which depends specifically on its mass. 


Figure 1. Simple evolutionary tracks for the 
internal conditions of stars of various 


masses 
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The general evolutionary trend of contraction, in- 
creasing the central density and temperature, is interrupted 
periodically by nuclear burning. The energy-generation 
history of a star is a succession of gravitational con- 
tractions which raise the central temperature of the star 
sufficiently to initiate thermonuclear reactions; the 
thermonuclear reactions transform a given type of fuel 
nuclei into heavier nuclei and release energy; the supply 
the given fuel nuclei becomes exhausted and the core 
resumes its gravitational contraction. The order of theme 
nuclear reactions is determined by the nuclei present and 
their charges. The larger the nuclear charge,' the higher 
its Coulomb barrier and the higher the kinetic energy 
(temperature) of the bombarding particles must be to 
penetrate .the barrier apd initiate nuclear reactions. 

A schematic sketch of ; th. energy history of a star is shown 

“ PigUre 2 ‘ Duri "« «»ol-r burning the temperature is 
almost constant. During gravitational contraction the iso- 
topic composition does not change. 

The most abundant element i s hydrogen, which also 
has the lowest charge, one. It is transformed into He 4 , 
releasing 6 x 10 18 -erg/g at temperatures above 10? * (Reeves, 
1965). Helium is transformed into C 12 at temperatures above 
about 10 -K and at slightly higher temperatures the carbon 
reacts again with helium to form o“. The amounts ofj carbon ’ 
and oxygen produced in the core during helium burninj depend 
on the central temperature and therefore on the mass if the 
star. The C and 0 curves in Figure 2 are the lower aid upper 


Figure 2. Energy history of a star (schematic diagram). 
Nuclear burning stages and the resulting composition 
of the core of the star are shown. Where two curves 
are drawn they represent the lower and upper limits 
of the range of nuclei produced. (H. Reeves: Stellar 
Energy Sources . Goddard Institute for Space Studies, 
NASA, 1963). 


Equilibrium Process 



H Burning 
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| 

I 

limits respectively. Carbon reacts with itself at temp- 

8 

eratures above about 7 x 10 °X; carbon burning produces 

nuclei in the range 0^ to Mg^ . The two curves again 

are the upper and lower limits. Neon photodisintegrates 

and oxygen reacts with itself at still higher temperatures, 
9 

about 1.4 x 10 °Kr- Neon burning predominantly produces 

"0^ and Mg^. Oxygen produces isotopes in the mass range 

28 

A = 25 - 32 with a strong peak at Si . The two curves 
show the lower and upper limits. 

The full chain of thermonuclear reactions does not 
occur in all stars. For a star of given mass there is a 
maximum central temperature attainable in a nondegenerate 
core. The exclusion principle requires that the average 
separation of particles be greater than the electron wave- 
length , 



,/2 m e 'kT 


, ( 2 . 10 ) 


where r is the size of a cube containing one proton and 
X. e = a/P and P = J2 m e kT. Using expressions (1.13) and 
(2.5) for p - and T, we must have 


1 > 


m 


'—'1/3 
ft VP 

1/3 /2" mTk? 


.0914 U 


-1/2 /MgA 


M 


V6 A, 



Thus the condition for nondegeneracy requires 


( 


a * is 


1/3 


( 2 . 11 ) 


The necessary central temperature for hydrogen burning is 
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T c £ 10 °K, 


so that mass and radius must satisfy the condition, from 
(1.13) and (2.5), 

7 

/*i jt \ / n _ \ T A * 

- 2.16 . (2.12) 


/M \ /R©\ 10 

^ Iw©)Uj 


4.61 x 10 ^ 


Combining these two requirements, (2-11) and (2.12), the 
minimum mass a star can have and burn hydrogen is 

3/2 m 


£L * 0.05 . 


(2.13) 


8 


For heliwn burning the central temperature must be 10 “K. 

The maximum central temperature occurs when the hydrogen 
burning shell has burnt its way almost to the surface, so 
W e can treat the core as a homogeneous star. The min im u m 
mass for helium burning is thus 


3/2 M 0 M a 

l c % * * 278 ° r % * * 18 * 


(2.14) 


The necessary central temperature for carbon burning is 

O 

about T = 7 x 10° *K. The minimum mass for carbon burning 
c 


is 


u 3/2 

c a*© 


(2.15) 


The necessary central temperature for neon and oxygen burning 

1 

is T ** 1.3 x 10^ °K. The minimum mass for neon and oxygen 
c 


burning is 


,3/2 -2£- i i o 

M 0 * i-’ • 


(2.16) 


Oxygen and neon burning are the end point of thermo- 
nuclear burning stages. Nuclear reactions among larger 

mass nuclei (further photodisintegrations and recombinations) 

9 

would occur in the temperature range of 2 - 4 x 10 °K'. 

However, at these temperatures the rate of energy dissipation 
by neutrinos (which are produced in the core and escape 
directly from the star) is so large that further nuclear 
reactions are unable to halt, but can merely slow down, 
the gravitational contraction. These reactions can, how- 
ever, produce nuclei all the way up to Fe^ , and the temp- 
erature is high enough to produce statistical equilibrium 

among the various nuclei. 


III. EARLY STAGES OF EVOLUTION - HOMOGENEOUS STARS 


Hydrostatic equilibrium and overall energy conservation 
determine the evolution of central stellar conditions. 

For more of the details of evolution, including the star's 
radius and luminosity, the mode of energy transport from 
the interior to the surface must also be considered. 

The equations of stellar structure — mass conservation, 
hydrostatic equilibrium, energy conservation, and energy 
transport — form a system of nonlinear differential equations 
which must be integrated numerically. It is possible, 
however, to obtain crude analytic stellar models by separat- 
ing the condition of hydrostatic equilibrium from the energy 
'transport . In the previous section, the condition of over- 
all hydrostatic equilibrium was expressed by the virial 
theorem. Now, since a more detailed stellar model is 
desired, we' assume an 1 analytic density distribution, 
namely, that the density in a star varies linearly from 
the center to the surface. (Cameron, 1963 ) . It is then 
possible to integrate the equations of mass conservation, 
hydrostatic equilibrium and energy generation through the 
star. Hence, together with the equation of state of an 
ideal gas, the run of density, mass, pressure, temperature, 
and luminosity through the star are determined. Also, 
the central density, pressure and temperature, and the 
total rate of energy generation are determined as a function 
of the star's mass and radius. Finally, the different 
modes of energy transport — radiative transport with Kramer's 
or electron scattering opacity and convective transport — are 
considered. The energy transport equation can be satisfied 
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at only one typical point of the star because of the 
approximation made in assuming a given density distribution. 
This gives a mass- luminosity- radius relation which gives 
the evolutionary track of the star in the Hertzsprung- Russell 
diagram. 

To summarize: Hydrostatic equilibrium and energy 

conservation determine the changes in the central stellar 
conditions, while hydrostatic equilibrium and the mode of 
energy transport determine the changes in the surface con- 
ditions — the track in the Hertz sprung- Russell diagram. 

Linear Stellar Model i 


Assume the density in a star varies linearly from 
the center to the surface. 


p.(r) = p c (l - -) 


(3.1) 


where R is the radius of the star. We call this a linear 
star model. The equations of hydrostatic equilibrium and 
energy generation can now be integrated but the energy 
transport equation can only be satisfied at one point in 
the star. The mass distribution is (from equation (1.4)) 


M(r) = J 4TTr 2 p(r)dr 


o 

T p c r ” 


4TT (1 -i|). 


4 R 


(3.2) 


Hence 

Thus 


M(R) = ^ ttP R 3 . 

3 c 


= 

ttRJ 


= 5. 6 j(JL\/*aY 


g/ cm' 


(3.3) 
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The pressure is obtained from the equation of hydrostatic 
equilibrium (l.l) 

P = P ^ P r GM(r)p (r)dr 

c Jo r 2 ' | 

where P c is the pressure at the center. Hence I 

p - p c - f a ?y ci - 1 £ + 1 i 

i 

i 

Applying the boundary condition P(R) = 0, we get j 


P = Gp 2 R 2 , 

c 3 6 c 


(3.4) 


hence 


24r 2 2.8 r 3 

R3 " R4 * 


* - 3* GpgR 2 (5 - ^5 

- ^-44 * 10 1S (k) 2 feVd 


2 3 

- 4.8 jz + 5.6 - 1.8 -gj) , (3-5) 


Assume that the radiation pressure is negligible; the 
temperature is then given by the perfect gas law equation (1.2), 

_ LlP 

T = — ■ 

kN Q p * 

where N 0 is Avagadro' f s number, the number of nucleons per 
gram. 


rr 




R' 


r2 

R 2 


* 3 6 kN 0 "c- 

= 9.62 X 10 6 ^(f^Cl - i - 3.S 

(3.6) 

We now know how the density, temperature, and pressure 


r 3 

+ 1.8 „ 

R3 
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vary throughout the interior of this linear star model. We 
have found the condition of hydrostatic equilibrium. The 
run of pressure, temperature, and density through the star 
is shown in Figure 3» We must now consider the condition 
of energy conservation. The rate of thermonuclear energy 
generation can be expressed in the form (Reeves, 1965) 


& - e G p 


k ( T \ n 

> [Jjrj ergs/g-sec # 

The luminosity of the star varies as (equation 1.6) 


= 1* 4TTr' 2 p(r') £dr* 

• J o 


For a linear density distribution the total energy gen- 
'eration is 

.R 


L =* f 4TTp(r)&r 2 dr 

« n 


SO 


l = 4TtR 3 e o Pc 


0 


16 


P ( 

tr & o' 


JL GuH l_^n M 


12 * T 


( 1+k+n ) 
n ( 3 k+n ) 


Thus 


L_ 

*<D 


, c cs , T ^-962^" „ f MN n+k * 1 /Ro'i n+3k 

35-58 e o I n u [tc] ^ 


(3.7) 


where 


I n - $1 x 2 (l-x)" +k+1 (1 + 2x - 1.8x 2 ) n dx 


has values cf the order of 10”^ or 10” ^ i and- T Q £^ is the 
temperature in units of 10^ °K. The energy generation and 
luminosity in a 1 ^ star are shown in Figure 4. 


i 



log />(g/cm 3 ) log T(°f<) log P(dynos/cm 



Figure 4 


Central Hydrogen Burning 






- 21 - 


Radiative Energy Transport 


Finally, consider the equation that governs the flow 
of energy through the star. First consider radiative energy 
transport. The appropriate equation for radiative energy 
transport- is equation (1.8) , 


T - „ 2 4ac T dT 

. L 4-TTr —7T- — — j — , 

• r o K p dr 


so the temperature gradient necessary to drive the radia- 
tive flux through the star is 


dT = _ 3*.. o L r 

dr 4ac T^ 4rr.r 2 


We consider two types of opacity (Cox, 1965 ) : (l) 

Kramer’s opacity^ 

x = h g P t ^ ) (3.8) 

which is a good approximation at intermediate internal 
temperatures, and ( 2 ) electron scattering opacity, 

x - * e - -20 (1 + X), (3.9) 

(Where X is the mass fraction of hydrogen), which is 

dominant at high internal temperatures. We also assume 

for convenience that all the energy is generated at the 

center of the star, so that 

L = L = constant, 
r 

For Kramer’s opacity 

dT = 3xa P 2 L 
dr 4ac 1^5 417 r 2 


(3.10) 
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Corapare this expression for dT/dr with the radial deriva- 
tive of T from the linear model. 


dT 

dr 


TT Q-tP c 

36 kN o 


R (5-38^ + 27 

R R z 


( 3 . 11 ) 


For our analytic model, these two expressions for the 
temperature gradient cannot be equal throughout the star, 
We determine the luminosity by equating the above two 
expressions at r = 0 . 5 R. 


, .2 4ac. T?/o /dT'V 

= - 4tt(t> , (— 

" 1/2 3k 0 Pi/2 \ dr 'l/2 ) 


where r^y 2 ~ 0 * 5 R > 


T l /2 


31 

T* 

288 kN 


Gup_ 2 
£ R 


( 


’ 1 / 2 , 

: dr\ 

, dr /l /2 


- 0.5 P, 


- 29rr Gu 

144 kN 0 p c R< 


Now 


hence 


3M 

ttR J 


_ 3 29/3 1 \ 6 * 5 ac /GHV* J 7.5 M 

L =TT 81 W ^ W W 


.7 . 5 


5.5 

R‘5 


J (3.12) 


where a = 7*57 x 10 ^ ergs/ cm^ deg^, c = 3 x 10 ^ cm/sec, 
G = 6.67 x 10“® dynes cm^/g 2 , k = I .38 x 10 ^ ergs/deg, 

“ 1.67 x 10”"^4 g t Hence, for population I stars 


ar<d H 


L 


- 0.5 


(3.13a) 
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5.5, 


7.5 / M V’7 R \ 

— = 28.6 \x — l j — I 

% V%>) W) 


and for population II stars 
L 


. _ 7*5 /M^ 5 /R^ 0#5 

— ■ 147 U • (tH 

i© y%>; vr© j 


(3.13b) 


Solar matter is ~ 2/3 hydrogen ' and l/3 

helium by weight. The mean molecular weight for 12 nucleons, 
of which 8 are hydrogen atoms and 1 is a helium atom';, is 

8 x 1 + 1 x 4 ! 


U = 


2X121 3 3 


number of particles 8 x 2 + 1 x 3 


12 

19 


- 0.632 


7.5 

H = 0 . 0320 . 


Hence 


L = 3.36 x 10 

For solar mass and radius 


33 


5-5 x-0.5 


ik\ ' (a 


er 


g/sec 


33 

Lq = 3-36 x 10 ergs/ sec. 


compared with the observed luminosity of the sun which is 

3 3 

L© ~ 3.89 x 10 ergs/sec. Thus the linear star model 

gives a result which is within 20 # of the observed value. 

The luminosity increases rapidly with the mass of the star 
and increases slightly with decreasing radius. 

When electron scattering is the dominant opacity the 
temperature gradient needed to transport the energy flux L is 

dT _ 3x a p L 

dr 4 ac T^ 4 rrr 2 


(3.14) 
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Determining the luminosity by equating this temperature 
gradient with the expression for dT/dr obtained in the 
linear model (3.11) at the midpoint r = 0.5R, gives 


which is 


t - .2 iac 

L - - 3> . e 

?!/;, f dT N | 

Pl/2 l dr A/2 ’ 


, = 29 2 ( 31 \ 3 

T ^2SS/> 

W •£ 

V k / H e 

(3.15) 

electron scattering 

dominates. 


L _ 178 

4 /M_\ 3 

(3.16) 

T "I ' -\r '4 

L o i T A 

( %/ ‘ 


The luminosity is independent of the radius and increases 
with mass, although less sensitively than for Kramer 1 s 
opacity. 

Equations (3*13) and (3*16) are the radiative mass- 
luminosity- radius relations for Kramer's and electron 
scattering opacity. The effective surface temperature is 
defined bv 


Flux = a T 


4 

eff 


or 


= 5*7 6 x 10 

Convective Energy Transport 


T e = (L/ 4 TTaR 2 )^ A A 


(3.17) 


The convective energy flux is (equation 1.10) 


f c = p c p w e , 

where w is the velocity and 0 the temperature fluctuation 

i 

of the moving convective element. Convection is extremely 
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efficient, and therefore all the energy to be transported 

can be moved by convection with only negligible adjustment 

in the superadiabatic gradient - \ — >. The 

Id r Vir/ ad ) 

energy flux is thus determined by the boundary layer of 
the convective region (Spiegel, 1965) . 

If the star has a substantial region with radiative 
transport, that region will determine the energy flux. If, 
however, the star is completely convective, the boundary 
layer determining the flux is the thin radiative photosphere 
surrounding the convective zone, where the energy must be 
transported by radiation since the material is becoming 
optically thin. The luminosity of the star is then. de- 
termined by the temperature of the gas at the point from 
which photons can escape from the star, F ra( j = crT^, so 


• 9 A 

L = 4nR“crx 


e > 


where T e is the effective surface temperature of the star. 

The depth in the star from which photons can escape 
nearly coincides with the transition point between the 
radiative and convective regions and occurs at an optical 
depth of about 2/3. The radiative temperature gradient 
drops rapidly as the density decreases, so- the temperature 
is practically constant from this point outward. We thus 
assume an isothermal photosphere and take the effective 
temperature as the temperature at the transition point 
between the convective and radiative regions (Hoyle and 
Schwarzschild, 1955 and Hayashi, Hoshi and Sugimoto, 1962). 
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We assume an opacity law of the fora 


x “ 


Then, since the bottom of the photosphere is at an 
optical depth 2/3, 

f * pdr fr Vr , 


and from equation (l.l) 


2 _ b 1 

3 *o T e „■ 


1 dP 

P g dr 3 

f p%p - 

J r (a— 1) g 


where P p ^ is the pressure at che bottom of the photosphere. 
.Thus, one relation between the temperature and pressure 
(or density) at the bottom of the photosphere is 


T p ', = f (a+1) 


( 3 . 18 ) 


This relation is the boundary condition for the star 

P ~ V “ | (a+1) ^ SS T - T aff • . 

This condition is just that the photon mean free path 
(xP) equals the scale height P/p g at the boundary so that 
the radiation can escape from the star at the effective 
temperature . 

A second condition on T and P .. can be obtained from 

e pn 

the condition for the boundary of the convective zone, namely. 


F C = f r 
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In the expression for the convective flux (1.10), let us 
approximate the velocity w by half the sound velocity > 


c “ 7 YkT/uH 


since c is an upper limit to the velocity. Also let us 
approximate pC 0 by y times the internal energy 

Sr 


U = 1 


kT 


UH 


- 4 p . 


Then the convective flux is 


? c = 5 pC w6 » $ Y | U 


(y — 

8 v ^ K / 


\ 'u { 

- v P T~ 


(3.19) 


The radiative flux is 


' F R = ° T e 


(3.20) 


Thus, equating (3.19) and (3.20), the transition point is 
given by j. 

3-5 

s T e • (3.21) 



The conditions (3.18>) and (3.21) can be combined to 


determine the effective temperature, which is 

1+3 - ' ', <1+a)/2 - T 2 ,-m 01 l/<W-3.Sft+«» 

te>J 

(3.22) 


= 0 


(l + a) 


;l+a 




u 


In the outer layers of stars the opacity is due primarily 
to H and is an increasing function of pressure and tem- 
perature, so a, b > 0. The H opacity is very temperature 
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sensitive, so b is large. Thus ^ e ff is nearly constant; 
it increases slightly with increasing mass and decreases 
slightly with increasing radius. 

The approximate power law form for the opacity ob- 
tained from the detailed opacity calculations in the ; 
region about 3500 °K is: ! 
For population I stars (X = 0.6, Y = 0.38,. Z = 0.02) 


x. = 6.9 x 10 


-26 T,0.7 m 5*3 


and for population II stars (X = 0.9, Y = 0.099, Z = 0.001) 

*- 6.1 X 10- 40 p°- 6 r 9 - 4 , 

where X, Y, Z are the mass fractions of hydrogen, helium 
and all the heavier elements respectively. The luminosity 
is found by inverting equation ( 3 . 17), 


= 4TTR 2 aT^/lQ = 


T e V (R_\ 

5.76 x 10 3 J 


.2 


(3.23) 


Then the effective temperature and luminosity, are: 
For population I 



and for population II 



-he effective temperature is less sensitive to the radius 
for population II than for population I stars because the 


- 29 - 


opacity is more sensitive to temperature. In population II 
stars, there are fewer metals with low ionization potentials 
to provide electrons to form H . The electrons must now come 
partly from the ionization of hydrogen which has a high 
ionization potential, so the electron pressure will be very 
temperature sensitive. 

In stars with high surface density, the relation (3.2l) 
between the pressure and temperature at the bottom of the 
photosphere is not valid, because in deriving it from the 
boundary condition Fq = we evaluated the convective flux 
by assuming that the temperature fluctuation is of the order 
of magnitude of the temperature itself. This assumption 
is valid only in stars where convection is inefficient near 
the surface due to low density and large radiative losses 
from the convective elements. In stars with high surface 
density, convection is very efficient and the temperature 
gradient in ohe convective region is nearly adiabatic 
throughout. In this case, the temperature fluctuations 
are much smaller than the order of magnitude of the temper- 
ature itself. 

For stars with high surface density therefore we go 
to the opposite extreme from the low surface density case and 
assume the temperature gradient is adiabatic throughout the 
convective zone. We may then use the adiabatic relation 
between pressure and temperature. In the interior. 


= K T Y/(Y ' 1) 


K T 


2.5 


P 


t 
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since (neglecting radiation pressure) y = 5 / 3 , except in 
the hydrogen ionization zone* For a fully convective star 


K — const . = P /t 
c' c 


2*5 


From the linear model (3*5) and ( 3 . 6 ), 


o 

GM“ 


P c 4rr R 4 


thus 


K 


T = — 

c 12 


GyH M 
k R * 


5 M -°* 5 R - 1 * 5 


= _5_/l2\ 2 °/R_\“^ -i, 

4rrl 5/ V'^ H / 

In particular, the above relation holds at the bottom of 
the hydrogen ionization zone. 

If we neglect the effect of hydrogen ionization, 
which reduces y , then at the boundary between the convective 
zone and the photosphere 

P ph “ K T e’ S 

with the same K as for the interior. This relation, 
combined with the optical depth condition equation ( 3 .I 8 ), 
gives the effective temperature 


T e = 


GM - ( 1+a )1 + 2.5(l+a)] 

•? (1+a) 75? K J 

2-5 


= fl 

(3 *0 


t H^j 2 ' 5 1+a ^ 2 . 5 + 1 . 5a 


fin . i 

Is V 12 ,’ ^ 

2.5(l+a) 1.5+0. 5a 1.5a-0*5) VCb+2. 5(l+a) ] 
xn M R 


I 
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For population I 

T 


2.6 x 10 3 |i° 


.445 ALt 194 ^ 0 ' 0576 

VMo) (.%> j 


_L 

L© 


— = O.O4I U 


!.78 f m _^°*776 ^ 2.23 


%/ V^o 


For population II 

T 


-3.01X103 (^f^- 0 * 98 . 

k - “•«” ^ 192 (it) 


0.686 


2.119 


The hydrogen ionization can, however, be treated 

exactly and we can relate K =. P . /T at the top of 

e P n 

^ <y r 

the hydrogen ionization zone to K = P^/T ^ 2 * 9 = P /T *' -at 
its bottom. Since the temperature varies adiabatically 
through the ionization zone, the entropy is constant 
across it. The effect of the ionization zone is to decrease 


d In T = ? - 1 

din P r 


so that the temperature will decrease less than the pressure 


going outward through the ionization zone. Then K q < K 
and will be increased.. The entropy per unit mass is 


c - 29 s F/ n , - % 5 x_ 

s - “TT ? < l+x+6 )- T kT 


H |_V 2 

+ x In 


In 


'_2trH') 3 / 2 


h V 


+ 6 In ( 


8 tth ^ 

h 2 / 




+ ( l+x+6 ) In 


( kT ) 5 / 2 d+x + 6 ) 


where y is the ionization energy of hydrogen, 6 = Y/4X, and 
x _s the fraction of hydrogen ionized. Evaluating S = con- 
stant above and below the hydrogen ionization zone, that is, 
for x = 0 and x = 1, respectively, gives 
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.= K = (1+6) 


r„ \ 3 / 2 /kX 5 / 2 " 1 ” 1+ 6 / *. \ (2 + 6)/( 1+ 6) 

■'“’to 0 J w 

. cm) . 


Thus the effective temperature is 


— l/ ( 1+6 ) r 

= £2 fl +a x 1 2 -49 /. 2 + _ 6 

ef f 5- 1 ' 7 "rP i 1+6 (l-« * TO" 2 


Cii (3 K o%' L 1 T 6 \1 

2 . 5 ( 1+ a r M \ 1+ °- S ( 1+ a) ( 2+6 5 /( 1+6 5 

x u (^m 3 1 

/r \1.5(l+a)f^ - 2 -) VCb + 2.5(1+.)] 

X : x \ 

\*o} 


V 1 • 53 x 10‘ 


2+6 _ . 
1+6 1 1_ra 


(3.27) 


-Again, in the high surface density as in the low surface 
density case, .the effective temperature is very insensitive 
to mass and radius. 


For population I 


, A .3 0.829 / M °* 27 ( M °* 288 

3.06 x 10 u (v ) VrJ 


= 0 


_. 0 3 . 32 / mn 1 * 08 / R _ y* 15 
6 U 1%) ^Ko) 


For population II 


3 0.59 JM \ °- 1 ^ 2 5 ( R n.0.204 

p W W ! 


*e - 3-75 x 10 4 ^ 

L 2 36/M\ 0 -”CS'l 2 ' 816 

k - oa785 u (*) U ) 


0 . 28 ) 


( 3 . 29 ) 


S umma r-y : The central 

linear stellar model are 


conditions of a star in the 


Pc = 5 - 64 &)(?) ’ 

P c - 4.44 X 10 15 , 


(3.3*) 


(3. S') 
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T - 9.62 x 10° n 

V 


ft) ({>) • 


( 3 . 6 ’) 


The total rate of energy generation is 


_L 

L, 


0 /. 9 62 \ n ' n /M \ n+k+1 f W >+3k 

- 3S ' S8 e o I n (t o(7) ) 4 (%j U/ 


where the rate of energy generation per gram is 


e - e o 0’ 


/it 


n > 

l o/ 


7 

and T /_•> ‘is in units of 10 degrees, and 
o(7 ) 

I = f 1 x 2 (l-x) n " rK ‘ 1 (l-2x-1.8x 2 ) n dx . 

n -i o 

The evolutionary tracks of stars in the Kertzsprung-Russell 
diagram depend on* the mode of ■ energy transport, which 
determines the mass-luminosity-radius relation. For fully 
convective stars, the luminosity is determined by the 
surface condition. Since the opacity is very temperature 
sensitive, the effective temperature is nearly constant, 
independent of the radius, and the track in the Hertzsprung— 
Russell diagram is a nearly vertical line. For stars with 
radiative energy transport, the luminosity is nearly independent 
of the radius and the track in the Hertzsprung-Russell 
diagram nearly a horizontal line. The changes in the stellar 
radius depend on the sources of energy and the internal 
structure of the star. 


A. PRE-MAIN SEQUENCE CONTRACTION PHASE 


The linear stellar model is now applied to. the pre— 
* 

main sequence contraction stage of evolution. A star is 
formed from a condensation of the interstellar gas that 
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is dense enough to become opaque to its own radiation. 

Then as the gas contracts its temperature will rise. As 
the temperature rises, the gas, composed predominantly of 
hydrogen and helium, is ionized. Much energy is necessary 
to ionize the gas, which means that the temperature cannot 
rise much above 10“ ‘K until the hydrogen is ionized. The 
ionization of the hydrogen and helium leads to gravitational 
instability, since the energy released by the contraction 
does not increase the kinetic energy per particle (the 
temperature) but goes into the ionization energy of the atoms 
Hence, the contraction of the gas does not raise the pressure 
sufficiently to permit the gas to remain in hydrostatic 
equiliorium, the ratio of specific heats y falls below 4 / 3 , 
and the collapse must continue. 

A stable star is not formed until the hydrogen and 
helium are almost completely ionized throughout most of 
the gas fragment. In such a contracting star, with the 
internal temperature of the order of 10 5 -' *K, ’ the opacity 
is so high that the radiative transport of energy i s impeded. 
Further, the extensive ionization zones increase. the specific 
heat and reduce Y to 'less than 4/3 throughout large regions 
of the star. Thus the adiabatic gradient 


f 


dT 

dr/ ad 


Y - 1 UH 
Y k . S 


will be small and the star will be unstable to convection 
throughout most of its interior. Its luminosity will then 
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be determined by the surface conditions. For a given star, 
the rate of contraction is limited by the rate at which 
energy can be radiated away 

1 dR _ LR 

R. dt a GM 2 ” * 

where a « 1. Thus the stellar structure with the highest 

luminosity will be stable and the condition for a fully 

convective star is that the convective luminosity exceed 

the radiative luminosity. Initially, stars that are not' 

too massive will be fully convective. Very massive stars 

(M > 12 Mq for population I and M ~ 16 >1 q for population IX) 

have > L cqrv and never pass through a fully convective 

-stage. Inclusion of the radiation pressure will, however, 

modify this result by increasing the convective instability. 

We first determine the conditions for the contracting 

star to become stable. The condition for stability is that 

all the hydrogen be ionized, 

kT > — 13*6 volts 
3 • 

or 

T > 5*2 x 10 4 °K 

The mean temperature of an homogeneous star is 



Thus the maximum radius of a stable star is 



= no 


f m \ 
* 1 % ; 


( 3.3 0 ) 


nrv» „ 
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effective temperature and luminosity of such marginally 


stable stars, as given by the fully convective linear model 
for low surface density ( 3 * 24 ) and ( 3 . 25 )* are 


e 

L 

T_ = 


h> 


, rt 3 — 0 . 2 OO 

3.15 x 10 u 


/mJ 

l M©, 


-0.089 


, 3 0.988 /m \ 1*644 

1 .07 x 10 u — \ , 

V ^0/ ' 

, _ ft 3 -0.187 /mV 0 * 0667 ") 
3.19 x 10 u 1 — 1 1 


Pop. I, 


1.152 x 10 3 (|) 


1.13 [_ 1,733 POP ‘ 


II, 


These relations give the starting point for the evolution 
of stars. 

As a star contracts, when fully convective, the effective 
temperature is nearly constant. The track in the H-R 
diagram follows the mass- luminosity- radius relation for a 
fully convective star (equations 3 •24*. 3 *25,. 3 . 28 ,^ 
and 3-29)1* For population I (X - 0.6, Y = o’.38, Z = 0.02) 



log 



- 7.236 log T e 
10.94 log T e - 


+ log 
1.874 



log 


r 0.843 log (1 + 28.345 

(low surface density), 

1 M \ (3 .31) 

(.%)- S ' 75 lo e “ * 39.77 

(high surface density) . 


These tracks are shown in Figure 5* For population II 
(X = 0.9, Y = 0.099, Z = 0.001) 


log 


log 


S3- 


11 log T e + 


;n 




0.821 log u + 41.485 

(low surface density), 

(X \ (3.32) 

- 13.8 log T e - 1.887 log^-J- 5.78 log u - 50.07 

(high surface density) . 


Figure 5* Hertzsprung-Russell diagram of the pre-main 
sequence contraction evolutionary tracks and 
initial main sequence. The tracks are labeled 
with the type of energy transport determining 
the direction of that portion of the track. 

Dotted curve is observed main sequence (Hayashi, 
Hoshi, and Sugimoto, 1962, and Schwa rzschild, 1957). 
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These tracks are similar to, but slightly steeper than those 
for population I stars. 1 

As a star contracts, its central temperature increases 

> . 

according to (equation 3*6). 


UH GM 


0 • l 

The increasing temperature increases the emission of 


radia- 


is in 


tion and reduces the opacity. A central core which 
radiative equilibrium will develop. When about half the 
star is in radiative equilibrium, the star will leave the 

fully convective path. The luminosity will now be deter- 

< i 

\ 

mined by the radiative flux, which is proportional to j 
(equation 1.8) i / 


1 dT 
K dr 


3 \ 

since T /p is approximately constant. The opacity decreases as the 

temperature rises. Thus as the star contraots, the lum- 
inosity will increase slightly, and must rise. The 

star will then move to the left in the H-R diagram. 

The radiative mass-luminosity- radius relation, for 
Kramer's opacity is given by equation (3*13) and the path 
in the H-R diagram will be 


logi 


/L \ fM \ 

( l ^) 83 °*8 lo £ T e + 1 os [ mq ) + ** log ^ “ 2*02 * 


If the central temperature becomes very high and the central 
density is low, the dominant opacity is due to electron 
scattering. Then the mass- luminosity- radius relation is 
given by equation ( 3 . 16), 


* 


I 


» 
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L . 178 . 4 

TfX w 


/M'? 

/ * 


A <>4- / vr \3^ / n \“j 

2.14 x 10 4 (1+X) 4 *u f|M ^ ) . 


This is the path followed by the massive stars. Typical 
radiative tracks in the H— R diagram are shown in Figure 5 


Time scale of Contraction 


The luminosity of a star is the rate of change of 


total energy 


AE 

At 


1 

2 At 


The gravitational energy is (from equation 1.19) 


- n 


G M z 
R 


so 




L pa 2 


GM* 

RAt 


Thus the time scale of the contraction phase is 

,2 


At - 

LR 


(3.33) 


- 1.59 x 10 7 years. 

Pre-main sequence ‘contraction times are listed in Table 1. 

B . CENTRAL HYDROGEN BURNING 

As a star contracts, its central temperature rises 
until it is high enough for hydrogen thermonuclear reactions 
to produce the energy radiated away from the star. At this 

point, the contraction stops arid the star spends most of 


i 

[i , 


Table 1 - Evolutionary Time Scales (Years) 


1 

Mass 
(Mo) j 

Population 

Pre-main 

sequence 

contraction 

Central 

hydrogen 

burning 

Hydrogen 

shell 

burning 

Central 

helium 

burning 

,7 

I 

2 X IQ 8 

„ ,„10 
4 X 10 

1 X IQ 7 

6 X IQ 7 

V 

II 

2 x 10 8 

5 x 10 

2 x 10 7 

4 x 10 7 

1 

I 

4 x 10 7 

9 

8 x 10 

6 x 10 6 

3 x 10 7 


II 

4 x 10 7 

9 x 10 9 

/' 9 x 10 8 

2 x 10 7 

2 

I 

3 x ‘ 10 6 

5 x 10 8 

2 x 10 6 

1 x 10 7 


II 

3 x 10 6 

7 x 10 8 

' 3 x 10 6 

9 x 10 6 

5 

I 

3 x 10 8 

3 x 10 7 

3 x 10 5 

2 x 10 7 


II 

1 x 10 6 

8 x 10 7 

4 x 10 5 

2 x 10 7 

7 

i 

l I 

2 x 10 5 

1 x 10 7 

2 x 10 5 

8 x 10 6 


II 

7 x 10 5 

4 x 10 7 

1 x 10 5 

8 x 10 6 

10 

I 

1 x 10 5 

7 x 10 6 

7 x 10 4 

3 x 10 6 


II 

, 3 x 10 5 

2 x 10 7 

4 x 10 4 

3 x 10 6 

15.6 

■ 

1 

6 x 10 4 

3 x 10 6 

2 x 10 4 

1 x 10 6 


II 

2 x 10 5 

1 x 10 7 

2 x 10 4 

1 x 10 6 


II 
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its lifetime burning hydrogen into helium. The locus of 
luminosity vs effective surface temperature of such stars 
(burning hydrogen in their cores and. still of nearly homo- _ 
geneous composition) defines the main sequence in the 
Hertzsprung-Russell diagram. 

The luminosity of a star is determined mainly by 
the thermal conductivity (radiative) of the stellar material. 
The central temperature is determined by the adjustment of 
the nuclear energy generation to maintain mechanical and 
thermal equilibrium throughout the star. Nuclear energy 
generation processes are very temperature-sensitive and 
thus nuclear energy sources play the role of thermostats. 

The radius of the star depends on the temperature and mass 
distribution. 

The basic features of the structure of homogeneous 
stars can be determined by dimensional analysis. The 
dependence of the central temperature and density on 
chemical composition, mass and radius is determined by 
the condition of hydrostatic equilibrium and the equation 
of state (from equations 1.15 and 1.13) 


T « nB U 
c R / 


M/R 3 . 


(3.34) 


The luminosity and radius are then determined by the energy 
balance. The equation for radiative energy transport is (1.8) 
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T , 2l6fl T J dl . 
L **. 4 TT r — - — — * 

3 Kp dr 


Assuming an opacity law of the form . 

a m b 


K o P ^ * 


the luminosity is 


-1 


L « * ( W 8) M 


4-b 3-a-b 3a+b 


(3.35) 


The rate of nuclear energy generation is (equation 1.6) 

L *» 4TrJ* epr 2 dr . 

Assuming the rate of nuclear qnergy generation per gram 
has the form 

e = e o p k T n , 

the total rate of energy generation is 

- / >n 1+k+n -3k-n 

L « e Q (uB) M R . (3.36) 

When the rate of energy generation equals the rate of 

energy loss (luminosity), then the dependence of the 

radius, luminosity and effective temperature, on the mass 

and chemical composition is 

, a . l/t , (n+b-4)/t (k+n+a+b-2)A 

R « (S 0 k 0 ) (HP) M , 

L oe K -( n+ 3k)/t e (3a+b)/t ^ 0 ^[n(4+3a) + 3k(4-b)]A 


x M 


[n(3+2a) + k(9-2b) + 3a + b]A 

( 3.375 


T 4 « K - (n+3k-2) A e (3a+b-2)A 

e o o 


x (uS )^ n ( 2+ 3a) + 3k(4-b) - 2b+8]A M Cn(l+2a)+k(7-2b)+a-b+4]A 

where t = n + 3k + 3a + b and b s 0 in the interior. The 
central temperature and density are 


i 

» 

5 

1 


T 

c 



(e 0 *o> 


-1/4 


<*o\>> 


-3A 
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, 0 ^(4+3k+3a)A 2(k+a+l)A 

y 

(3.38) 

(U8) “ 3(n+b “^)A M -2(n+b-3)A 


Thus the radius, luminosity, effective temperature and 
central temperature increase with mass, and the central 
density increases with mass for the p-p chain, n = 4, but 
decreases with mass for the CNO cycle, ti m 18. 

The main sequence is the locus of points in the 
luminosity-effective temperature diagram 



4 

4 

4 


n(3+2a) + k(9-2b) + 3a + b , m , 
n(l+2a) +' k(7-2b) + a- b + 4 log T e + const * 

i 0 g t + const. (Kramer’s) 

3n + 15.5 e 

(3.39) 

n ~ + log T e + const. (electron scattering). 


The central temperature of a contracting star is 


T c - 9.62 x 10 7 4 . 

Hydrogen burning starts at about T c = 8 x 10^ °K. 'Thus 
a star will start generating its energy by nuclear reactions 
when its radius is 

• 1 ' 2U (^) ' ’ (3 ' 40) 

Stars of small mass, M ^ 2^, burn hydrogen by the 

7 

p-p chain at a temperature around 1.5 x 10 °K. The rate 

of energy generation is approximately 

* ‘ 6 o p (l.l x io'V ^s/e-sec, 

2 
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Massive stars, M Z, 2Mq, burn hydrogen by the CNO cycle 

7 

at a temperature generation around 2 x 10 °K. The rate 

of energy generation is approximately 

18 


e = e c p 


{ttw) eres/e - 


sec j 

/ 

e 0 - 4Si x H x CN0 , 

The energy generation rates for the linear model are (from 
equation 3-7) 


= 4.98 x 10~ 3 u 4 


(ib) (r 3 ) p-p;chain^ 


i- = 1.157 x 1<T 4 u 14 Uj-'f 6 [j?! 7 CNO cycle . 

The properties of stars on the main sequence— —burning 
hydrogen in their cores — are: For the p-p chain and 

Kramer r s opacity 

-0.538/m \°- 076 9 

0.312 U ! TJ — \ , 


R 
^0 

T“ = 49.1 






»’■” (a 


5.46 


T = 2.18 x 10 4 u' 

e 


2.21 ( M X 1 * 058 


© 


ftV 


log(-r— )= 5-16 log T - 0.74 log U - 20.7 , 


7 1 o / m\®*923 

3.05 X 10 7 u 1 * 54 \ 


l.QS( m \°*769 

1 ‘ U>) 


p c = 186 u- 

For the. CNO cycle and Kramer’s opacity 


(3.41) 
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. ' 0.395 / M \ 0 *697 

°- 451 14 he) 


43 - s u?-3 (^) 


o 
5.18 
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*© 

L 

T e - 2.34 x 10 4 u 3 (^ 


log 


/JLl 

W 


T <c “ 


Pc = 


5.948 log T g - 2.39 log \X - 24.36, 

l .« x ^ u “^(|) 0,364 , 

8 u-0*455 .MX ' 0 - 909 (3.42) 

6S.Su 


Stars switch over from the p-p chain to the CNO cycle 

7 

at a central temperature about 2 x 10 °K, which occurs 
at a mass of about M “ 2M©. For CNO cycle and electron 
. scattering opacity 

\ 0.765 


II 

0.454 U 0 - 588 1 

L 

% 

• 4 /M \3 

112 u ( y M Q j 

- T e " 

2.77 x 10 4 H 

ios iky 
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e 
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2.12 x 10 7 U 
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' „ -1.765 

60.3 U 


3 /m \°*368 

V%>) 

.76 log U - 3 

i f *£\’ 235 

K) ’ 

. -1.294 


(3.43) 


Stars switch over from Kramer's to electron scattering as 
the dominant opacity for mass M > 3 ^ for population I 
and M > 2 ^ for population II. 

The evolutionary tracks for different mass stars 

I 

are shown in Figures 5 through 8 . During the pre-main 


Figure 6. Evolutionary tracks of population I stars in 
H-R diagram during pre-main sequence contraction. 
The main-sequence is also shown. .Dotted curve is 
observed main sequence. 
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Figure 8. Evolutionary track in H-R diagram of star at 
one solar mass. Solid and dotted curves are from 
analytic models. Dashed-dot curve is results of 
calculations using Henyey method by D. Ezer and 
A.G.W. Cameron: This conference, p. 
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sequence contraction the stars contract to release grav- 
itational potential energy to supply the radiative energy 

losses from the surface of the star. The radius of the 

* 

star decreases. The direction of the track is determined 
by the mode of energy transport: convection with low surface 

density, convection with high surface density, radiation 
with electron scattering opacity, or radiation with Kramer's 
opacity. Stars are fully convective when they first become 
stable, except for very massive stars M > 12 (population I) 
and M > 17 M© (population II) . Stars become radiative when 
the radiative luminosity is greater than the fully con- 
vective luminosity. 

The main-sequence is the region of the H-R diagram 

where central hydrogen burning occurs. Here the central 

temperature is high enough for the hydrogen thermonuclear 

reactions to supply the energy radiated away. There are 

three sections of the main sequence with different slopes, 

depending on the mode of energy generation and the type 

of opacity. Because the linear model is not sufficiently 

centrally condensed, the main sequence is shifted to lower 

effective temperature and slightly higher luminosity than 

0 

obtained from accurate calculations. The radius must 
shrink in order to raise the central temperature to high 
enough values to generate the luminosity. 

j 

i 

( 

Convective Core 


A star which is generating energy at its center by 
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a very temperature sensitive thermonuclear reaction (all 
processes except the equilibrium p-p chain) will have 
a convective core* The energy generation region is very 
small, so the luminosity increases very rapidly with radius. 

The flux F *» L/4rrr will then be extremely large, since 

! 

the radius is very small, which forces the radiative 

t 

temperature gradient to become superadiabatic in order 
to carry the flux. This causes instability to convection. 

f 

The boundary condition for the convective core is 


fdT\ J dT\ 
Vdr/rad \dr/ad 


(3.44) 


(*£) = 
\dr/ad 


1 1 dP _ 1 T GM r p ! 

(N+l) ad P d? ‘ (N+1 ^ad T r 2 ' (3 * 45) 


and 


(N+l) 


ad 


32 - 243 - 38 
8-60 


(3i46) 


where 


( 


r dT\ 

drjrad 


B - Pg/P 


3 HQ L _ 

16a t3 4Ttr z (N+l) rad P dr 


1 T dP . 

p j_ > (3 *47 ) 


where 


where 


(N+l) 


rad 


16ttcG(1-8)M t 
x L 


(3*48) 


(1-9) = P rad /P aT 4 /P . 


Thus the condition for convective instability is 

(N+1) rad S (1,+1 >ad • 


(3.49) 
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Expressed in another form, the boundary of the convective 
core will be at 

= 3 HP L r 

q M ' ^'ad lonacG M 

1 nL r 


(N+l) 


ad 16 ttcG (l-B)M 


(3.50) 


For a convective core to exist, the effective poi., tropic 
index must be N ft{ j and decreasing inward at some point in the 
star (Naur and Osterbrock, 1953), i.e., at the core boundary 


d In (X+l) rad 


* 0 . 


d In r 


a — b 

Assuming H = x p T , then 


d In (N+l) rad = d In T d In _ d In P _ d In h _ d In L r 

d In r d In r d In r d In r d In r d In r 


where 


" P£r - (1+a) ] 


d In P + d In M r 
d In r d In. r 


d In L r 
d In r 


- - (^7=? - V + B - W, 


N+l 


U = d In M r _ :4Trr^P 


V = - 


W = 


d In r 
d In P 
d In r 

d In L r 
d In r 


M r 

GM^p 

rP 


Expand M r , P, L r , and T about their central values, 

M r " | "P c p3 * 

2 2 2 

P = P c _ | ttG p c r , 


T = T„ - 


1 T 


L_ ■ 4 TT P- 

r 3 c c * 
1 2 


c " N+l p AP - T c " U+l 2 1X0 TT Pc 


2 2 
Pc '** • 


At' the center, U c = 3 > V c = 0 and W c 88 3, so ■ 
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d In (N+l)^^ / d In r a 0 at the center. Thus the condition 

4 +b+a 


rad 

for a convective core 

D . d In (N+l) rad 
d In r 




N+l 


- a 


- 1^ V + u - w 


* 0 


becomes 


dD s 0, 

dV 


since D c = 0 and V increases outward. 


Evaluate dU/dV and dW/dV at the center. U “ 3 t so dU = 0. 
Thus we must develop p and M to higher order. 


. _ UH P _ , N 2 r PcZ 2 \ 

p " H T ' p c (1 “ N+l 3 P~ r ) 

c 


I ''C'- N+l 3 

- Ocd - wTT c r2 >. 


2 Pc 

where C = — ttG — — . Then 

3 


M r = 4 n ,f^pr 2 dr = ^ p c r 3 (1 - | Cr 2 ). 


Then 


u - 3(1 - C r 2 ). 


so 


dU 

dV 


V - 2Cr 

3 N 


5 N+l * 


Now consider W. 


L r " 4 tt&o Jo p 


1+d -V 2 


T r dr 


d v 

assuming an energy generation rate of the form t = £_p T . 


0 = - ill C r2) and T = T o (1 * ill C r2) - 
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Thus 


- ™ 0 °™Tc So (1 - ^ Cr 2 ) (1 - * Cr’)r*dr 

Cr 2 ) 


4TT a . ^1+d m V _3 „ 3 y ‘~KN(l+d) „ 2, 


&o t; r° (1 - j 


T" ~~ K c *c * - -fe N + 1 


4 TT 3 , 3 V + N ( 1 +d ) 2 \ 

— e c (>o rJ f 1 - j ~+ i Cr '• 


Then 


[in (1 - a v cr*)l 


so 


and 


Thus 


W “ 3 + d In r LJ -“ VA " 5 

_ , 6 y .±_N(l + dl « 2 

3 “ 5 N + 1 » 


... _ 12 v + N(d+1) r . 

dW - 5 N + 1 Cr dr ^ 


dV = 4Cr dr 


dW 

dV 


_3 y + N(l+d) 
5 N + 1 


The criterion for the existence of a convective core 
is thus 

S “ TDOTT (3v + 3Nd + 5Na + 5N - 5b - 15) * 0, 

(3.51) 


where 


N + 1 = 


32 - 243 - 33 
8-68 


e - e «V , 


H •- K Q 0 a T" b 


For Kramer's opacity and 8—1, d = 1, this becomes 


v s 4*3 
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For a gravitationally contracting core & ~ T; for electron 
scattering opacity the condition for a convective core is 
then N > 2.4 which occurs for 8 = 0.75. For 8 = 0, no 
temperature or density dependence of & is needed 
in order to have a convective core. 

Assuming the existence of a convective core and L — L 

V 

at the core boundary, its size is given by 


9l 


(Ntl) a( j ^ L 
1-8 16 TTcGM 




(3.52) 


The size of the convective core depends on the mass of the 
star only through the radiation pressure. Consider the 
special case of negligible radiation pressure. 


1 = 3P = 3 k 

1-3 a 


3k 


P c 


T 4 aiiH ** aUH "r? i 

c 


3 

since p/T is approximately constant through a star. Then 


“ < N+1 >ad * 


3kL 


T J l6tracGuHM 
c 


Apply dimensional analysis to this expression. 



For electron scattering x. = and L ~ 8^ M^, so 

q x ~ 3 (N+l) ad . 

For Kramer's opacity K = k q d T j ' ; 



T ^ g -2/(n+2.5) M 6/(n+2.5) 

L ° _ g(7n+22.5)/(n+2.5) M (5n+15.5)/(n+2.5) 


q l 


~ 8 


1 + 6/(n+2.5) 


Now consider the case when radiation, pressure is important 


8 < l. Electron scattering will then be dominant, so; 




_ (N+Dad L 


and 


i— ~ — TT — 3 

1-8 T • 


1-3 

-4 w”2 


M 


M 


L a 4 ..2 
j while m ~ ° ” 


Thus for massive stars 

~ (X+Dad - 

which increases by a factor of 2 as 8 decreases from 1 
to ($. Thus in all cases depends on the mass only 

through the radiation pressure. 

In the case of negligible radiation pressure,, the 
size of the convective core can be found explicitly by 

using the linear model 

15 _}£_ i© ( *c°c 

q l 128tt GHct M© IWvM / UT c 3 * 

For electron scattering 

= . 22 (l+X) , 

and for Kramer’s opacity 

q, = 1 . 6 5 ( 1+X ) Z . 


IV . ADVANCED STAGES of EVOLUTI ON - INHOMOGENEOUS STARS 

A star spends most of its life burning hydrogen into 
helium in its core. The advanced stages of evolutxon 
comprise the star's life after central hydrogen burning. 
When the hydrogen in the core is completely transformed 
into helium, the core of the star contracts and heats up. 
The rising temperature enables hydrogen thermonuclear 
reactions to occur in a hydrogen burning shell source 
surrounding the core. A star in this stage is composed 
of a helium core, a hydrogen burning shell source and a 
hydrogen envelope. Depending on its mass, a star may 
proceed on to helium, carbon, neon, and oxygen burning. 

If the star is massive enougfy the core continues to con 
tract and heat up, until, at about 10 8 °K, helium burnxng 

thermonuclear reactions occur in the core. 

As a star evolves, each nuclear burning process starts 
first in the core and burns outward as the star heats up. 
Thus, a star that has passed through several nuclear 
burning stages will be composed of concentric shells of 


the products of the different processes, with a hydrogen 
envelope on the outside and a core of the products of the 
last nuclear burning stage through which the star has 
passed. Figure 9 illustrates the shell structure of a 
star that has passed through all the nuclear burning stages 


Figure 9» Schematic shell structure of a massive star at 
the end of nuclear burning. The star is assumed to 
have passed through all the nuclear burning stages plus 
tv approaching equilibrium among the nuclei in the core. 
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We first consider some general properties of stars in 

advanced stages of evolution. The evolution of stars is 

» 

towards greater central condensation. Stars contract and 
increase their central density and temperature. This con- 
traction is occasionally interrupted (but the evolutionary 
trend is not altered) by nuclear burning in the core of 
the star.. 

The increasing central density as a star evolves, to- 
gether with the existence of nuclear burning, shell sources, 
causes the development of large radii and extended envelopes. 
The large radii are caused by increasing central condensation 
that is, increasing central density but decreasing envelope 
density. The degree of central condensation is measured by 

U = In M(r) _ 4TTr P _ p (r) 

d In r M(r) p r * ' 

where o r is the mean density interior to r. 'Since 

d In r = — dlnq , 

U 


the radius is 


In 


1 R = 7r d ln + ln R i * (4.1) 

qi u i 


where and refer to the core-envelope interface. Now, 
from equation (1.15), 

Gu c H M X 


T £3 — -- ■ - 

c k R^ * 


SO 


is , 


(4.2) 


? 

where M, is the mass of the core. Thus the stellar radius is 


In R - 


In 


(4.3) 
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+ j: 


u 


d In q + In 




. l c/ q x 

The larger the central condensation, the smaller the U near 
the shell source and the larger the stellar radius. All 
stars in advanced stages of evolution have extended en- 
velopes . 

We digress now to discuss the nondimensional variables 
U, V, and N+l: 

U 


_ d In M(r) = 4nr 3 p 


d In r M(r) 

d In P _ GM(r )o = 3 


d In r 


rP 


P (r) 

- > 

Pr 

GM(r)/r 

2 - i p/p 


(4.4) 


xr+T 2 d In P = 16 Ttac GM(r) T 4 

d In T 3 P* ^(r) 

At the center of a star U 3 9 V-* 0, and at the surface 
U 0, V- 4 ®* The polytropic index N varies between 1*5 
for a convective region and infinity for an isothermal 
region* Also 


d_ 

In 

_T = 

V 

d 

In 

r 

N+l 

d 

In 

_P _ 

NY 

d 

In 

r 

N+l 


(4.5) 


Thus the r — dependence of the physical variables is given 
in terms of U, V, N+l by 


M(r) ~ r^ , 

P ~ r-V 
T ~ r-VAN+l)^ 
p ~ r - N V/(N+l). 


(4.6) 


From hydrostatic and thermal equilibrium, the physical 
variables r, M(r), P, and T must be continuous throughout a 
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star. A discontinuity in pressure would entail an infinite 

acceleration, and a discontinuity in temperature would 

* 

entail an infinite energy flux. At a composition dis- 
continuity then, the density will be discontinuous, but 
p/u will be continuous. Thus the ' continuity conditions 
on U, V, N+l are 

— , Z xL (N+l) Continuous. ( 4 . 7 ) 

U U r 

The dependence of the radius on the central condensation 
U can now be evaluated approximately. 


In R = .fq d In q + In , ( 4 . 1 1 ) 

where is the radius of the base of the envelope. The 

integral may be evaluated approximately by expanding U 

about its value Uj at the base of the envelope (from equations ( 4 ^ 6 ) 


and (1.4) ) 


so 


and 




NV/(N+1) 


= n /, NV Ar ^ 
P 1 U “ N+l RT >> 


4TT 


U = 


M(r) = + 4'n'r^P2Ar , 

(Kl 3 + 3Rl 2 4r) Pi (l - £r) 


M 1 + ^R-^p-jAr 


4 *rr Rj p 1 


(1 + 3 - 


NV 

N+l 


Ar 


- “V Rj 


u, (4 - - u,) — , 

1 N+l l 7 R-, * 

Ar _ AM 


(4.8) 


R 1 4ttR 1 ^P 


As , 

• U 1 
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thus 


AU 


(3 - 


NV 

N+l 


- U x ) Aq 


(4.9) 


The main contributions to the integral for the radius, 

(4.1), come from those regions < q < where U but 

not Aq is small; that is, not near q s 1, Let U = + ^q-Qj) 

where a = 3 - NV/(N+l) - U . Then 


A In R » J* q ° 


q i u i + a (q-qi) 


I dq 


( 4 . 10 ) 


- ( 1 \ m ( 2 * "i ( V 

V u i ■ aq l/ l/ 1 ! U 1 + 


For « aq^ 


For » aq 1 


For U » Ctq^ 


A In R « 


A In R m 


A In R 


aq x 

1 

aq x 

1_ 

U-, 


In 


\ U ] 


qp - qj 
<Jo . 


m [ 2 *\ . 

^1) 


Thus for great central condensation, small U^, the radius R 
is large. 

i 

The increasing central condensation in advanced stages 
of stellar evolution is caused by the increasing central 
density in conjunction with the existence of a shell energy 
source. Increasing the central density increases the pressure 
gradient dP/dr ** -pg. However, the core luminosity is less 
than the total luminosity, the core tends toward an iso- 
thermal condition, and the temperature varies by less than 


T'« 0 


1/3 


Thus the density gradient in the core increases. 
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* 30 i/p decreases, and the stellar radius increases. 

The composition discontinuity between the hydrogen envelope 

* 

and the helium core causes a decrease in and so by 
a factor of U c /w e , and also contributes to increasing the 
stellar radius. 

Although stellar radii tend to increase during the ; 
advanced stages of evolution, their actual magnitude de- 
pends on the detailed structure of the star. There is 
a general empirical rule for determining the variation of 
a star’s radius: The direction of expansion or contraction 

in a star is reversed at every nuclear burning shell source 
and unaffected by any inactive shell. The reversal of 
expansion or contraction of a nuclear burning shell source 
is due to the thermostatic nature of a nuclear energy 
source. A star adjusts itself to maintain a constant 
temperature in the nuclear energy source, which causes the 
radii of the nuclear burning shell sources to tend to 
remain nearly constant. The mechanism is similar to 
that which keeps a main sequence star in equilibrium. 
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Xf the radii of the shell sources remain constant, 
the contraction of a zone between two shells, for instance, 
means that the density at the inner shell of the zone in- 
creases but that the density at the outer shell must de- 
crease, since the mass and volume of the zone remains con- 
stant. Thus the density at the inner shell of the next 
outer zone is decreasing and that zone is expanding (see 
Figure 10) . 

Consider the zone between two shells of radii R Q < R^. 


R 


- r 


R-, 


-o * ^ 1 

Let m be the mass of this zone and assume R, » R . 

i. o 

mean density of the zone is 


M X 


- M, 


^ (R][ - R l) 


3m 

» >T 

4TT Rj* 


The 


(4.H) 


Thus 


ARj 

R 1 


1 

3 


Ap 

P 


( 4 . 12 ) 


We also assume the radiation pressure is negligible so 8 « 1. 

Consider what happens when the radius of the inner shell 
changes. Suppose R Q changes by AR 0 . If the shell at R Q is 
not nuclear burning, its properties vary in a manner that pre- 
serves hydrostatic equilibrium, that is approximately homo- 
logously. Then, by equations (1.13) and (1.15) > 


T <x — 
R 


I? 


Figure 10 
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so 


&T 0 


ARq 

T 


R 

0 


O 

Ap 0 

= - 3 

AR o 

A> 


R o 

**o 

" - 4 

ARq 

P o 


R o 


I 

, (4.13) 


Then 


ARj = 1 Ap 0 _ AR q 

*1 ' 3 p 0 Ro ' 

Thus when the inner shell is not nuclear burning the outer 
shell’s radius changes in the same way as the inner shell’s 
radius, and the shell has no effect on the expansion or 
contraction. 

If the shell is nuclear burning the structure of the 
shell initially changes homologously. However, due to 
the change in the rate- of energy generation, there • is an 
additional, no nhomologous, change in the structure. The 
change in the rate of energy generation is 


e N = e o (p+Ap)(T+AT) n = (1+ ^- + n -^) 

= «N 0 (1 - (n+3) , 

so that 

where 1^2 \ ± s the initial change in R 0 . Initially, this 

v R o J± 

net change in energy is deposited (or removed) where it is 
generated and the material heats up (or cools down) . The 
temperature changes until the fractional change in luminosity 
(rate of removal of energy from the region) is equal to the 


A&N 


ft 


- (n+3) 


N 


r AR c 
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fractional change in the rate of energy generation. Since 



gives the additional, nonhomologous change in T Q . This 

additional temperature change produces an additional 

* 

pressure change (besides that produced by the initial 
homologous transformation), 

n+3 / AM 

7-5 Uo Jl 

An increase in pressure produced by a contraction of the 
shell will push the shell back out; a decrease in pressure 
produced by an expansion of the shell will allow the shell 
to fall back in. The pressure must return to its equilibrium 
homologous value -and the shell must move back in the direction 
from which it came' according to the homologous relation 



AR 

R 


i AP 
4 P 


The secondary correction to the radius of the shell is 
hence 




n+3 /ARp\ 
4x7.5 V R 0 /i 
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The secondary correction to the radius is thus 



( 4 . 14 ) 


There is, therefore, a strong restoring force on the radii 
of nuclear" burning shells, tending to keep them constant. 

The changes in density, if the radii are precisely con- 
stant, can be found from the linear model. 


»(r)=p- (p 0 - Pl ) 


R 1 - R o 


where R ft is the radius of the inner shell and R. the radius 


of the outer shell. Then 


4TT 


* - M- - f- V CPo - - RJ]. 

R 1 -R 0 


‘1 ' “o ~ 3 "1 Lw ° ■ rT^I V 4^1 "O' 


The change in the mean density is zero. If R Q « R^, then 


A ?1 


Ac, 


(4.15) 


and changes in the opposite direction to p 0 . Therefore, 
sin che radii of the nuclear burning shells tend to re- 
aa... constant, the sign of the change in the density will 
alternate from one shell to the next. Thus the direction 
of expansion or contraction is reversed at a nuclear burning 


a n e -i_ 1 . 
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Apply the general rule for stellar radii changes 
to the various stages of evolution. During the pre-main 
sequence contraction stage, the core is contracting; there 

0 

are no shells, so the whole star is contracting. During 
the hydrogen exhaustion phase, the core is contracting; 
there are no shells, so the whole 'star is contracting. 

Then a hydrogeq burning shell is ignited, the helium 
core continues to contract, but now there is one shell, 
so the envelope expands. When the central helium burn- : 
ing commences the core expands; there is one shell, so 
the envelope contracts. These structural changes are 

illustrated in Figure 11. The structural changes during 

> 

the stage of helium burning are illustrated in C. Hayashi: 
"Advanced Stages of Evolution, " this conference, p. j . 

A. CENTRAL HYDROGEN DEPLETION 

We now consider in some detail the evolution of stars 
from the depletion of hydrogen in the core to the onWet ' 
of helium burning in the core. 

The depletion of a nuclear fuel in the core of a ; star 

. . ' , • i 

and the ignition of a shell source is a process which j 
'changes the basic structure of a star. We can there- I 
fore not construct an analytic model for this phase but 
only give some of its general properties. 

During central hydrogen burning, the luminosity 
of a star increases due to the increasing mean molecular 



Figure 11. Schematic diagram of the changes in stellar 
structure from the pre-main sequence contraction! to 
the onset of central helium burning. ' ! 
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weight as hydrogen is depleted in the core. Assuming 
the homology relations for homogeneous stars are still 
valid, in- small mass stars where Kramer's opacity is 
dominant from (3.13) and (3 .35) 

L « ( U 0) 7 - 5 f ' 

while in massive stars where electron scattering is < 
dominant from (3.16) and (3.35) ‘ ! 

L « (u0) 4 . 

The mean molecular weight increases by about a factor of 
2 as hydrogen is consumed. 

The energy generation rate has the form 

. n 

e - e o X H X 2 * (J^) , 

where X 2 is X H for the p-p chain and is Xcuo for the CNO 
cycle. The temperature exponent is n « 4 for the p-P 
chain and n w -17 for the CNO cycle. As the hydrogen 
concentration in the core decreases, the central temp- 
erature must rise in order to maintain the rate of energy 
generation. The radius of the star will therefore tend to 
shrink (see equation 4 . 2 ), 

. R « M/T 
• c 

The tendency of the radius to decrease due to the 

increasing central temperature is counteracted by the 

tendency of the radius to increase due to the growing 

\Xn + 

composition inhomogeneity Which decreases ^1- 

at the bottom of the envelope*. 


-63- 

The p— p chain is less sensitive to temperature and 
more sensitive to hydrogen poncentration than the CNO 
cycle. The central temperature wilL thus increase much 
more in small-mass than in large-mass stars. During • 

central hydrogen burning in small-mass stars, the rapidly 

j 

increasing central temperature nearly balances the grow- 
ing composition inhomogeneity and the radius stays nearly 
constant. In massive stars, the central temperature 
rises only slightly and the radius increases due to the 
composition inhomogeneity. The evolutionary track of 
a star in the H-R diagram during central hydrogen burning 


is towards higher luminosity. For low-mass stars, where 
the radius is approximately constant, the track’ is nearly 
parallel to the main sequence. For massive stars, where 
the radius increases, the track turns off the main sequence 
to lower effective temperatures. 

The equation for the consumption of nuclear fuel is 


dX 

dt 

dX 

dt 


e 

E 



radiative zone ; 
convective zone^ 


where X is the concentration of fuel nuclei and E is the 
energy released per gram of fuel consumed. This equation 
can be solved for the time scale of central nuclear burning 


M c 

At « y 2 E AX > (4.17) 

where L/M c « 6, the mean rate of energy generation, E is 
the energy release per gram, and AX « 1. Lifetimes of stars 
near the main sequence are given in Table 1* 
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B. HYDROGEN SHELL BURNING j 

\ 

As hydrogen is exhausted in the core of a star, the j 

central temperature increases in order to maintain the j 

rate of energy generation. The temperature farther out in J 

the star is then increased and the rate of hydrogen burning I 

It 

outside the core (where the hydrogen has not been ex-. I 

hausted) is therefore increased. Thus a shell burning source j 

is ignited. 

When hydrogen becomes nearly exhausted in small mass 
stars generating energy by the p-p chain, the central temper- 
ature has already increased and raised the temperature in 
surrounding regions of higher hydrogen concentration 
surf iciently to produce hydrogen thermonuclear reactions 
u..ere. When hydrogen becomes nearly exhausted in massive 
suars, the central temperature has not yet increased much 
due to the high temperature sensitivity of the CNO cycle. 

The energy requirements of the star must still be met by 
the core, so the central temperature must now increase 
greatly. This causes the radius of the star to contract 
and its track in the H-R diagram swings to higher effective 
temperatures. Eventually the decrease in X outruns the 

c 

n 

increase, in T and the rate of nuclear energy generation 
c 

in the core decreases. The core then starts to contract 
and release gravitational energy to supplement the de- 
creasing rate of central nuclear energy generation. The 


gravitational contraction raises the central temperature 
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1/3 


, and the shell temperature and ignites the 


shell source. The more massive the stax> the larger j the < 

size of the initial convective core and the farther 'out 
from the center lie the hydrogen rich regions. Ther^ the 
temperature in the hydrogen rich shell will be lower', the 
ignition of the shell source will be delayed, an<| the grav- 
itational energy release will supplant nuclear energy 
generation as the star’s primary energy source. Eventually 
the contraction will raise the temperature enough to ig- 
nite the shell source. Summarizing, as hydrogen is ex- 
hausted in the core of a star the temperature increase's, 

I 

nuclear energy generation in the core decreases, and a| 
hydrogen-burning shell source surrounding the core is 
ignited. 

When hydrogen is exhausted in the core and a shell 
burning source is set up, the pressure distribution in 
the core is initially similar to that of a homogeneous 
star. The value of U at the outside of the shell, U^+, 
is then decreased because: (a) The composition discontinuity 

U c /u e 2 reduces Oj+/p c and so by a factor of 2, and 

(b) when nuclear energy generation in the core ceases, 
the core tends to become isothermal. The reduced temper- 
ature gradient increases the density gradient, which re- 
duces Pi+/o c somewhat further. ^ 


The U - V locus of a star is given by 


l 


d In 

_U „ 

3 

- u - SSL , 

S 

1 

d In 

r 


N+l 

(4*18) 
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The points on the V - U curve with horizontal or vertical 
tangent are given by 

U + V/(N+l) -1=0 (horizontal), 

(4.19) 

U + NV/(N+l) - ,3 = 0 (vertical). 


These two lines intersect at the point 


U = 


N - 3 

N - 1 J 


V — 2 — — j * (4«20) 


Thus for N > 3 , the intersection point is in the physical 
region and there is a loop point corresponding to r -» • . 
Typical U — V curves for homogeneous and inhomogeneous 
stars are shown in Figure 12. 

For an isothermal core, N = «, so the U - V curve 
has a loop point of U = 1, V = 2 . The maximum V thus 
occurs for U = 1 and is somewhat larger than 2. Thus 
for an isothermal core 

U 1+ ea iis U « 0 . 5 • — -v 

U 

An isothermal core, if too large, however, cannot 
support the weight of the envelope. The critical size 
of an isothermal core can be found from the virial theorem 
(McCrea, 1957), 


Figure 12 


Centrally Condensed Envelopes 
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(v— i) — + .i £ 
VY X) V 3 V 


whore U and V are now the internal energy and volume. 
For an isothermal sphere the internal energy is, from 
equation ( 2 . 3 ), 


Y-l UH 


T M , 


and for a sphere of uniform density the gravitational 
energy is, from equation ( 1 . 19 ), 


3 GM^ 
5 R 


The pressure at the boundary of the isothermal core is 


therefore 


p = 3_ J£_ TM 

4TT (iH 'P 


3 1_ GM 

5 4ti F 


There is a maximum pressure consistent with:the equilibrium 


virial theorem, which is given by 


4TT nHR 


12 GM 
2 On -.R' 


Thus there is a critical core radius 


crit . 


= 4_ Gq r H ^ 


( 4 . 21 ) 


with stability possible only for R c0re 38 R cr it • 


maximum possible pressure is 


max i 6 n 


(1S\ 3 f rtA 4 1 

(^4 / V.UH j G 3 M x : 


, (4.22) 


which- decreases with increasing core mass. To determine 

the limiting mass of an isothermal core this must 

max 

be compared with the pressure necessary to support a star. 
For the linear model, (equation 3.5), 
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p = 5_ GM_ 

c 4TT R 4 


Also for the linear model, (equation 3*6), 


T l~ 


T = — 


5_ GuH M 
21 k R 


Thus the condition for a stable star, .that an isothermal 
nondegenerate core can support the surrounding envelope, is 


or 


P > P 

max ~ c 


£ 0.16 i 

M 


(4.23) 


Accurate calculations (Schonberg and Chandrasekhar, 1942). 
give < 0.1. 

If the mass of the core is below the isothermal core 
limiting mass (Schonberg-Chandrasekhar limit), the core 
becomes isothermal and the central temperature may decrease. 
In massive stars, the core exceeds the Schonberg-Chandrasekhar 
limit and gravitational contraction begins when nuclear 
energy generation ceases to support the star. In small 
mass stars^the core is initially below the limiting size, 
but shell burning adds material to the core until in this 
case, too, the core exceeds thfe Schonberg-Chandrasekhar 
limit . 

In all stars, therefore, to support the weight of 
the envelope the pressure gradient in the core must in- . 
crease. This raises the central density and greatly re- 
duces U 1+ = 3P 1+ /P C leading to very extended envelopes. 

The increased pressure gradient is achieved by two methods: 
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For small mass stars the electrons become degenerate 
and their degeneracy pressure greatly increases the pressure 
gradient. For large mass stars the core contracts rapidly, 
producing an increased density gradient and a nonzero 
temperature gradient, both of which combine to increase 
the pressure gradient . 

The envelopes of stars in advanced stages of evolution 
are therefore characterized by great extension, low density 
and small U near the shell at the base of the envelope. 

That is, the envelopes have a centrally condensed structure, 
with the density . increasing rapidly inward due to the 
large pressure gradient at the edge of the core, but the 
mass remaining nearly constant as r - r* shell from above, 
ihe greater the central condensation the larger the stellar 
radius. The most centrally condensed envelope structure is 
P as r -» 0. Since the mass, M r , must remain finite, 

0 06 r_a where a < 3. Then do/dr is finite, so from equation 
(4.5) V is finite. Also, since dT/dr ± 0 , from equation 
(4.5) N is finite. From equation ( 4 . 4 ), since M r is ap- 
proximately constant, V « (rT)” 1 , S o 

T 06 1/r • (4.24) 

We can now determine the limiting values of the nondimensional 
variables at the base of such an extremely centrally con-, 
densed envelope. Since r^g^ is very small, the values 
at the base of the envelope will not be too different from 
their values in the limit r - 0. For N < 3 , the limit as 
r - 0 is, from equations ( 4 . 6 ) and ( 4 . 24 ) 
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B-0, V-N+l. ( 4 . 25 ) 

Thus U ■* 0, the limit of extreme central condensation . 

The radial dependence of the physical variables is 

-(N+l) _ -1 


, T « r 
-N 


(4.26) 


P * r 

P oe r 

For N > 3» the limit as r -» 0 is a loop point given by 
equation (4*20). In this case U > 0 and the envelope 
is not so centrally condensed. The radial dependence 
of the physical variables is 

P ~ r -2<N+D/(N-l) ( T _ r -2/0f-l)' 

-2K/(N-1) (4.27) 

P ~ r . 


We now determine the effective polytropic index at 
the base of a centrally condensed envelope. In terms of 


A 


nondimensional variables 


V 


II 

GM 2 

4TTR 4 

> 

T = t 

UgH 

GM 


k 

R * 

cr 

II 

u 

X 

M , 


r — x 

R, 


equilibrium equations are 

d P = 
dx 

tx z 

>> 

CD 


^3. = 

dx t 


(4.28) 


St, 


v.-. . re t •' si/ u e , and the flux equations are 


v 


~ = _ c P 

dx K x 2 t 8 * 5 

. _ c p 


dx 


(Kramer's opacity), 

(4.29) 

(electron scattering). 
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For Kramer's opacity, combining equations ( 4 . 28 ) and ( 4 . 29 ), 

2 — 

dt®* s 4.25 % » 


so, near the surface and near the shell at the base of the 
envelope, where U is very small and • the mass fraction 
q is nearly constant, the polytropic index is 

N = 3.25 • 

Thus at the shell 


0 


-2.89 

r 


V* 

Similarly, for electron scattering. 


(Kramer's). ( 4 * 30 ) 



_q 

4 C E 


> 


so, near the surface and near the shell the polytropic 
index is 


N = 3- 

Thus the density distribution at the shell is 

D ~ r (electron scattering). 

(4.31) 

The only envelope model which can be readily solved 
analytically is p(r) ~ r ” 3 . This, as was just shown, 
corresponds to the limiting case of extreme central con- 
densation for both Kramer's and electron scattering opacity. 

The internal structure will be well represented by such a 

k 

model, bpt because it is too centrally condensed the stellar 
radii will be much too large. To calculate the radii a 
somewhat less centrally condensed model should be used. 
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Xnhomogeneous Analytic Stellar Mode] 


We now construct an analytic model of a star with 

\ 

one shell using a linear density distribution .in the core 

-3 

and an r density distribution in the envelope. 

1 • Core 

In the core assume a linear density distribution 


o(r) - o_ - (p. 


- R X • 


(4.32) 


where is the central density and p^ is the density at 
the shell. Then the mass distribution in the core is, ~ 
by equation ( 1 . 4 ), 

.r 2 

M(r) = J 4 rrp(r)r dr 


4TT 3 
r 


[o c - | ( 0 C - Pi) (4.33) 


and the mass of the core is 


V- 1 R 1 3(o c + 3D t)’ 


(4.34) 


This relation can be turned around to give the radius of 


s-(a.f(sr 


W* C 

°c + 3 o-jj 


-1/3 

(4.35) 


The pressure in the core is determined by hydrostatic 


equilibrium, equation (l.l). 


P (r) - F„ - af 


- V- f <*c 2r2 a - id- £=> id-^-^i 


M(r)p (r) dr 

h 1 - f" ,;+ 


I 
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At the boundary of the core 

p -JL 
1 




so 

p 


k tt 2 2 Pi Pi 

P l- T l + JS Gp c R 1 <* + 10 ^ + 9-f=a) . 


Pc 


Pc' 


Thus the pressure in the core is 


P(r) = p l- T l + 16 G °c V IS + 10 + 9 

u c 


Pi- 


' 24 £a + 28(1 - 77* £3 - - 2 77 + (4.36) 


Pi- A Pl-_ } r 

c M c ^1 


and the central pressure is 


P C u<j H P l- T l + 36 Gp c R 1 t 1 + 2 + T * 8 ~~2 }• (4.37) 

n 


k 

‘c 1 

c 

For a perfect gas, with negligible radiation pressure, 
the temperature is, by equation (1.2), 

T (r) = ^ . 

k p (r) 

Thus the temperature in the core is 


T(r) = [1 - (1 - |1=) f-]- 1 [ £I=.,t. + !i- 202 p R ' 

Pc R l P C 1 36 k p cM 


n Gu c H 

[j + 10 £1= +.9^ - 24-=^ + 28(1 -^=)£L 

Pc Pc 2 Rl* p c r7 

A 1 2 


(4.38) 


and the central temperature is 


c 


Pc" R 1 


u c Pi 

p; 


m , JSL p 2 n (. , .. Pl- 

T 1 ^ 36 k R 1 Pc* 1 + 2 — 


• ° 1 - 

+ 1.8 A ) 


(4.39) 

2 2 


Ue ^ T 1 + °-” * 107 u c(^) 2 »c(1 + 2 + 1.8 4 

VR©/ P c n| p 2 


i 
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However, when the core is degenerate it is assumed to be 

} 

isothermal, so T c ** in a degenerate core. j 


-3 


2. Envelope 

In the envelope assume an r ° density distribution 

3 

(4. 40) 


0 (r) 




Then the mass distribution is, from equation ( 1 . 4 ), 


M (r) * M + 4 ttp 1 R 1 3 f ~ j 

a. j. Rj 


V 


where is the mass inside the shell. Thus the mass dis- 
tribution in the envelope is 


M (r) = M 1 + 4tto 1 R 1 3 lnCr/R^, 

/ 

M = M x + 4 ttc 1 R 1 3 ln(R/R 1 f. 


(4. 41) 


The pressure is- determined by hydrostatic equilibrium, 
equation (l.l). 


P(r ) = P G f MlWll dr 
J R 1 r z 


= P 1 - Go^ 3 J* [M 1 + 4np 1 R 1 3 ln(r/R 1 )] dr. 


SO 


P(r) - 2 1 - jG Pl (1 - (£ l)\ - I Gp^ Rl : 


Th 


+ ttGp^Rj 6 -pf (i + In ^ ) . 


e boundary condition P(R) =0 determines Pj 


P(R) - 0 - P x - iG 0l ^ (1 - (j i)‘ 


) - -^GOi R x ' 


nGo i 2 + ln ' 
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so 

P 


i ~ iGOl ) + 4 Go i R ! - ttGPi 2 2l^(i + In 2-) . 

x \ / R7 | 


Thus the pressure in the envelope is 
P(r) 


<*.s[(?r- (?r] 

>.V [?;)■- (?r *.(?•)•>. ti ijl. 


(4.42) 


and the pressure at the shell is 




.2ft 4 

\R } ' 4 

4 


“ - (?r 


- ~.v 


) y 


(4.43 ) 


Note that the pressure is proportional to r“^ except near 
the surface. 

The temperature for a perfect gas with negligible 
radiation pressure is, from equation (1.2), 

T (p) = £ Zj£) •. 

k P (r) 

Thus the temperature in the envelope is 
GU„H 


T(r) 


4k 
+ trp 


1*1 {?- ? fty * - ? - ^ 

- - P (?Y i-X) ■ 


( 4 . 44 ) 


Xot-e that the temperature is proportional to except 

near the surface. The temperature of the shell is 


“ ft <>•(?) V <4 .-(MW 5,.}] 

.577 X 10 7 U e + .3066 X 10 7 U e(^) P 1 » 


( 4 . 45 ) 


The stellar radius is extremely sensitive to the 
degree of central condensation* for our envelope density 
distribution p ~ r"3, the radius is obtained from the 
mass relation ( 4 . 4 I), 



Thus the radiup depends exponentially on 0 c /t>±» This 
leads to extremely large radii, much larger than are 
observed. This is to be expected, since this envelope 
corresponds to the maximum degree of central condensation. 

Consider now less centrallj' condensed envelopes with 

\ 

density distributions 


P (r) 


fn \ n 
(?) 


(1.5 < n < 3 ) .. ( 4 . 47 ) 


The mass distribution 


in the envelope is then 


M (r) = M 1 + 4 np 1 R i n J'^ r 2-n dr 


= M + 


£L_ « D n r-3 ~ n „ 3-n 


1 3 


~ P ! R ! - V 


and 


M - *1 * ^ pa" (R 3 '" - R! 3 '") • 


Thus the radius is 


Rx + f (3-n)(M-M 1 ) ^ 1 /(3-n) 

3-n \ 4 tt pj R T n J 


R M - 


(4.48) 



log />((j/cfYi 3 ) log^E^K) log Ptetyfios/em 




Figure I4 
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degree of central condensation. As a rough approximation 
for all starsy we choose n “ 2\ 5 • Then i , 

R_ _ 2 *1 +/ «0 V f*V (X-o- ) 2 f R A~ 5 n ' 2 
V : % [W) Vife) (l<,l) (id' " 1 | 

-5 -2 >;i 


2 + 5.54 X io“ 2 (5 


is • iH 


The internal -structure of the inhomogeneous model ] 
is shown in Figure 13 • Figure 14 shows the tremendous 
degree of central condensation of the mass as compared 
with the homogeneous model. 

We now turn from the hydrostatics to the energy balance 
•in the envelope. The rate of thermonuclear energy gen- 
eration in the shell is, equation (1.6), 

T * (* R 2 h \ n 2 

. L " 4TT6 0 j p /--J r 2 dr, 

1 ' L o' • • * 


where we have assumed a nuclear energy generation rate per 


gram of the form 


e - e 0 o 


Then 


ItS- 


l - 4ne oPl 2 0l) d n • 

- ^vd$ n (>-( zD • 


iH£l_ e 0 . 

.n+3 0 


'Thus the rate of thermonuclear energy release from a shell 


source xs 


T - e 

n+3 




(ZT-HW- <*•»> 

The energy generation is confined to an extremely thin shell 
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source as shown in Figure 15 • 

The luminosity for radiative energy transport is 

2 l6ff T 3 ~ b dT' 

L = - 4 nr -*rx“ ~r~ » 

3 H 0 pl+a dr 

assuming an opacity law of the form 


* = K Q p a T b . 


The temperature gradient determined from hydrostatic 
equilibrium is i 


dT 

dr 


= - T 


R 1 

iT 2 


Thus the radiative luminosity is 

„ 64 tto „ m T 3_b 

L = — o Vi ^ • 

Evaluating the luminosity at the shell gives 


L = 


64 TTO 
3 *0 

64 TTCT 

rp 4 - -V 

n Tl 

1 Pl 1+a 
m 7-5 

R 1 * 

3 V 

R i 


4 

64 TTCT 

R. li- 

3 

1 »1 


(Kramer's), 

(4.51) 

(Electron Scattering), 


where for Kramer's opacity a *= 1, b — -3«5, and for electron 

scattering a = b = 0. Thus the radiative luminosity of 

the envelope is, for population I (X = 0.6, Y ** O. 38 , Z = 0.0 2 ), 

.3 


4.25 x 10 




(Kramer 1 s) 


U Q 


3 /r \ t .4 (4.52) 

2.18 x 10 !■— ) ^( . Z l . (Electron Scattering) 


Figure 15 
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and for population II (X— 0.9, Y— 0.099, Z = 0.001), 

2.27 x 10 4 (^Jl) (Kramer’s] 

(4.53) 


L 


0 


1.84 X 10‘ 


(*i) M 
1 %/ 0 2 


Lzl 


(Electron Scattering) # 


The effective temperature is given by equation ( 3 . 17 ). 

'.r 

For a fully convective envelope, the luminosity and 
effective temperature are determined by the surface condition, 
equations (3-24), (3.25), ( 3 . 28 ) and (3.29). The track is 
the same as for pre-main sequence fully convective con- 
traction, but traversed in the opposite direction. 

The time scale of evolution is determined by the rate 
of release of energy, 

t dE . , _ AE , . . 

L = dt s ° At T • (4.54) 

The time scale during stages of core contraction is de- 
termined by the gravitational energy release, 

2 


AE = i |Afl| = i 


GM 


c 


RR' 


AR 


The luminosity of the core is determined by the opacity 

(usually electron scattering since the temperature is high, 

~ 10 °K) and temperature gradient. In stars with degenerate 

* I 

cores, the central portion where degeneracy is strong is 

isothermal and the temperature drop occurs in the outer 

nondegenerate portion. The luminosity of a contracting 

core will therefore be assumed to be 

_ - 170 a. 4 /V 3 

I© ' 


^ - 179 ^(k) 


1 
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Thus for a contracting core 


At 




358 


/Me' 




3 ) 


N-l 


AR/R© 


(r/r©)(r7r 0 ) 


( 4 . 55 ) 


where R c is the previous and R c the current core radius. 
The amount of material added to the core by the hydrogen 
burning shell during this time is 

At L 


AM. 


so 


Aq x 


AM 

— 


e H X < 


% 




\ Mo L 

7M" X 




At , 


(4.56) 


Evolution During the Hvdrogen Shell Burning Phase: 


The- evolution during the hydrogen shell burning phase is •: 

’ ! | 1 

toward greater central density and temperature and larger stellar 
radii. The central density, o c , is chosen as the parameter 
labeling the course of evolution, since during the con- 
traction of the helium core, p.^ increases monotonically . 

A sequence of models, with increasing p c describes the 
course of evolution. It is necessary to choose an initial 
core size to start the sequence, since the details of the 
setting up of a shell source cannot be followed analytically. 

When p c » pj and R^ « R, the stellar structure can 
be expressed as an explicit function of q^, M, and o c . The 
core radius is, from equation (4*35), 


' 178 'h\ 1/3 (<>cS l/3 

• 178 \K) £ 5 *) 


R o 


04 . 57 ) 
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so shrinks with increasing p c . 

The central temperature for a nondegenerate core is 
found by substituting equation ( 4 .57) for in equation 
(4*39) and neglecting the first term. Which is negligible. 



(4.58) 


The shell temperature is found by substituting equation 
( 4 * 57 ) for in equation ( 4 . 45 ) and neglecting the second 
term, which is small^ 



For a nondegenerate core 


(4.59) 


V T i ■ i - 67 • 


Xn small mass stars, M < 3 - 4 Mq, the core is degenerate 
and isothermal, so 


T 

c 


*1/3 


Thus the core temperature increases as p c . The shell 
density is determined by the energy balance, luminosity 
= energy generation rate. The energy generation rate per 
gram is assumed to be of the form 


e - e Q pf- 


(0 


and all constants are evaluated for the CNO cycle at 

T o = 2 x 10 7 °K, so e o = 451 X H X CNQ and n - 18 . The -fcofcaJ 

energy generation rate is given by equation (4.50). 

The luminosity depends on the opacity and the energy 

transport mechanism, ~and two cases are considered: Rad- 

iative transfer with electron scattering opacity, equation 
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(4.51) , and convective transport with low surface density, 
equation (3*23) with equation (3*22)* For electron scat- 
tering 

' *V m C l*e 


-(n-4)/3 /MiVT (n “ 3) /^p c 

' UJ Uo3/ 1 


(4.60a) 


where 


. . f38. 

1 C42. 


2 (population I) 

8 (population II) 


I 


C 2^e 


For a convective envelope 

, P A-jl 1 f ( 2n+13 )A-4 o1 

(l- qi ) 2 LP=W. q * L 3A-8 J 

1-1 lU-uol («.60b) 


• */ 



_ 1 C (n-13 ) A+ 40 I 

Z L 3 A- 8 J * 


10 


where A *= b- + 2.5(a+l) for an H" opacity law of the form 


* - ^pV, 


and 


'7 


■’{ 


11 (population l) 

15 (population II) , 


C: 


5 (population I) 

6 (population II) 


The radius of the star is given by substituting equation 

( 4 . 57 ) for and .equation ( 4 . 60 ) for in equation ( 4 . 49 ) 

* 

and neglecting the small first term* For electron scattering 


R_ 

*0 


^ 


(kp 9)/9 


2 (4n-27)/9 2(n-4)/3 


x (lSv 


\ (?n+3)/9 


where 


[S: 


( 4 v 6 la) 

. | . 'i 

. 1 ? v. ,• ;; 

'j ; j • 1 ’ 

21 (population I) il iV-M 

17 (population II) * . < 
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For a convective .envelope 


r _ c „ PUSH 

*© 4 e W (4.61b) 

VlO^j 


d-q^ 1 


fl&l 




where 


t 


25 (population I) 


4 ( 7.0 (population II) 

Thus the stellar radius increases with increasing central 
density. The luminosity of the star is given by substituting 
equation ( 4 * 60 ) for and equation (4*59) for T^ in equation 
(4*50). For electron scattering 

L _ „ (n+ 8)/3 (m 1 \ <2 " +21)/ Vp c ^(" + 3)/9 

Wl 


- ,a t) &) 


(4.62a) 


m .(*1.1 x 10^ (population I) 

U .j 


where 

^ (. 8 .4 x 10 ^ (population II) 

For a convective envelope 

L 


1 ; “ Cft,le fcd 


(4.62b) 


x [tl- qi ) q 1 <n " 1)/3 ] 4 t C (n 3A-a A ' 4) l ( 


where 


p 2.7 

V.20.3 


(population I) 
(population II) 


The effective temperature of the star is found from equation 
( 3 . 17 ) with equations (4.62) and (4.61). For electron scattering 


where 


Ma 


o 

x U 
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-(2n-13)/l2 . .( 6n . 7S )/ 3 j 6 

(i-Ql) q x | 

-(n-8)./4 , 0 x -(n+l)/l2 , ^ 

s (-^\ ‘ , (4.63a) 


C: 


107 


24 x 10 4 (population I) 

55 x 10 4 (population II) 


For a convective envelope 

- w -faf&rl 


C 8 w e 


Ol) 

[,-HS'‘“ >/ ‘] foW (-&) 


(4.63b) 




where 


'8 


ft 


5.11 x 10^ (population I) 


62 ' x 10 J (population II) 


2 

For radiative envelopes, the luminosity increases, (~ 0 ), 

c 

and the effective temperature decreases, (~ *P C ), with 

increasing central density. For convective envelopes^, 

4*5 1 

the luminosity increases rapidly, (~ P c ), and the ef- 
fective temperature decreases slowly, (~ p c ^*^), witl^ 
increasing central density. The mode of energy transijort 
in the envelope switches from radiative to convective ! 


when the convective flux becomes larger than the radiative 
flux . 

The tip of the red giant sequence occurring in small 

mass stars is determined by the onset of helium burning. 

Helium burning commences in the center of a star when 

8 

T c « 10 *K. In stars of small mass with degenerate 
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cores the core is nearly isothermal. For an isothermal 
core 


T - T - 3.24 x 10' 

C JL 


1 . “• Gte) /3 (ifr) 


1/3 


(4.S9 1 ) 

so the central density at which helium burning commences is 

p -3 f MA" 2 

^7 - 29.5 U 4 

10 3 4 


(ST 


Thus the maximum luminosity at the tip of the red giant 
branch, where the envelope is convective, is, from equation 
(4.62b), 

,* „ -*1^1 


L 

io 


const .. (9.54) 

■cir” 







(4.64) 


At T - 9.5 x 10 ‘ 
o 


(i). 


max 


r 


22 x 10 


94 x 10- 


c / M \0.108 n-aA 1,189 :I 

W tql 1 ) (Population I) 




Thus the luminosity at the tip of the red giant branch at 
the onset of central helium burning is very insensitive to 
the mass of the star, and is about two orders of magnitude 
higher than obtained from accurate calculations (see Figure 20 ) 
due to the absence of any temperature gradient in the core. 

The evolutionary changes in the central conditions are 
shown in Figure. 16 . For low mass stars the core isj iso- 
thermal and the central temperature is constant or liay even 
decrease slightly when the hydrogen in the core is exhausted 
and a shell source is ignited. The central density | increases 
with nearly constant central temperature until the increasing 




Figure 16* Evolution of central conditions during pre-main 
sequence contraction,' central hydrogen burning, 
helium core contraction and oentral helium burning. 

The solid lines and shaded regions are from the analytic 
models, the dashed dot lines are interpolations* 

The dashed lines are from Hayashi, Hoshi, and Sugimoto, 
1962 . 
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s luminosity along the red giant branch causes the shell 

1/3 

temperature rises much more rapidly than p c until it 

8 ' i 

approaches 10 *K and helium thermonuclear reaction^ are ' 
ignited. The energy released by helium burning in the 

i 

degenerate core raises the central temperature, without 
affecting the density, until the material becomes nonde- 
generate. The core then expands reducing the central 
temperature and density. Stars with masses less than about 
3-4 ^ develop degenerate cores. \ 

In massive stars an isothermal condition does not 

| , 

develop. The core contraction provides an appreciably 
part of the star's luminosity from the beginning of j 

t 

• hydrogen shell burning. The central temperature and density 
increase, with T c increasing slightly less rapidly than 

1/3 - 

p c • 

The evolutionary tracks of stars in the H-R diagram 
during hydrogen shell burning are shown in Figures 17 to 
21. The stars' move to 'the right' in the H-R diagram be- 
cause their radii are increasing. The tracks depend in 
their grossest features on whether or not the star is. small 
enough to develop a degenerate core. Those stars that 
develop isothermal degenerate cores must evolve to much 
higher central densities and much greater central con- 
densation than those that do not. Thus very small mass 
stars develop very extensive envelopes, .which are there- 
fore fully convective and very luminous* These form the 
red giant branch (Figure 20). Intermediate mass'stars 


I 




s 







Figure 17. Evolutionary tracks of stars in H-R diagram during 
hydrogen shell burning with helium core contraction. The 
nature of the energy- transport mechanism in the envelope, 
which determines the slopes of the tracks is shown. The 

shaded area is the region where stars have just started to 
bum helium into carbon in their cores . : 


Hydrogen Shell Burning 
Population I 
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develop radii large enough to develop fully convective 

1 i 

•i 

but not such extensive envelopes) and their luminosity 

■i ' \ 

does not greatly increase (Figure 2l)». The very- massive 
stars do not develop a very great central condensation 

g 

. before their Central temperature has reached 10 °K, ; so 

they do not develop convective envelopes before helium 
burning. 

At the beginning of hydrogen shell burning the lumin- 
osity, except for very small mass population I stars, 

- x is much too low because we have taken the most centrally con- 

s 

densed model throughout and have not allowed the degree 
of central condensation of the envelope solutions to ; grad- 
ually increase. The temperature falls off extremely rap- 
idly outside the shell and the hydrogen shell burning region 
is therefore very thin, covering about 1 % instead of an 
initial 10 % of the mass, as found in accurate calculations. 

- The total amount of energy generated is therefore too small. 
Since the degree of central condensation and the thickness 
of the shell are constant, the luminosity in our models in- 
creases during evolution. Accurate calculations show, (how- 
ever', that the shells are originally much thicker than ours, 
and the narrowing of the shells, due to the steepening ’temp- 
erature gradient and the exhaustion of fuel on their inside, 
counteracts the rising shell temperature and the luminosity 
stays fairly constant. 

The. time scales for evolution during the hydrogen shell 
burning - contracting helium core stage are given in Table !• 



Figure 18 





log L/L© 






Figure 20. Evolutionary track in H-R diagram of 1.2 M© 

star. Solid lines are from analytic models for pre-main 
sequence contraction (PMSC), central hydrogen burning (H), 
hydrogen shell burning red giants (RG), and central 
helium burning (He) . Dashed line is from Hoyle and 
Schwarzschild, 1955 and Selberg and Schwarzschild. 
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Figure 21. Evolutionary track in H-R diagram of 7Mo star. 
Solid lines are from analytic model for pre-main ! 
sequence contraction (PMSC) , central hydrogen burning 
... (H), hydrogen shell burning (HSB), and central helium 
burning (He). The dashed curve is from Hoffmeister, 
Kippenhahn, and Weigert. 
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Summarizing: The cause of the extended envelopes 

of hydrogen shell burning stars is their central conden- 
sation; the cause of their central condensation is as 

P 

follows. As the hydrogen in the core is exhausted, the 
density and pressure distribution in the core change only 
slightly, but a chemical composition discontinuity de- 
velops and the density of the outside of. the core (at 
the shell) is decreased. Since U c /U e = 2 and p/u is 
continuous across the shell, the density at the outside 
the shell is halved. This decreases the ratio of the. 
density at the shell to the central density, i.e., increases 
the central condensation. When the core hydrogen is 
exhausted and the thermonuclear energy generation occurs 
in the shell, the core tends toward an isothermal state. 

Then the density gradient in the core increases and this 
further increases the central condensation. An isothermal 
nondegenerate core cannot, however, support, the weight of 
the envelope when the mass of the core is larger than 
about 10$ of the mass of the star. Then the pressure 
gradient in the core must increase in order to support the 
weight. This further increases the central condensation, 

P center/p shell* to very large values. 

C . CENTRAL HELIUM. BURNING 

When the central temperature of the helium core is 

i 

8 

raised to about 10 °K, helium will begin to burn at 

the center of the star. If the core was degenerate, a 
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v _ helium flash will occur because the pressure of degenerate 
matter depends only on the density, not the temperature, so 

that the energy released by the onset of helium burning will 

1 

increase the temperature without a corresponding increase 
in pressure. The increased temperature speeds up-tl^e 
reactions, which further increases the temperature, until 
the temperature is high enough for the matter to become 
nondegenerate. The rapid increase in the helium reaction 
rate continues until kT in the central degenerate region 
rises above the Fermi level and the perfect gas law again 
holds. In nondegenerate material, increasing the temper- 
ature increases the pressure, which causes the core to ex- 
• pand, thereby reducing the density and temperature and 
damping the reaction. The core will then settle down to 
burning helium at a much lower density and slightly higher 
temperature than at the onset of the flash. In stars with 
nondegenerate cores, there is no flash; the process of 
adjustment is 'small and occurs smoothly. 

A star burning helium at its center will be much 
more centrally condensed than a main sequence hydrogen 
burning star . The density at the shell where the com- 
position discontinuity, and possibly hydrogen burning, 
occurs is much less than the central density. Thus the 
core may be treated as a separate star with the density 
but not the temperature, going to zero at its surface. The 
luminosity of the core is determined by the balance between 
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the radiative energy transport and the helium energy 
generation rate. This energy balance determines the 
central temperature, which together with hydrostatic 
equilibrium determines the central density. The radius 
of the core is determined by the density distribution,, 
which we assume is linear. Thus for the model of a ' 
star burning helium at its center, assume a linear density 
distribution in its core, a p ~ r“3 density distribution 
in its envelope and treat the core as a separate star. 

The helium energy generation rate is 

e 3<* - C 12 - (4.66) 

for 

Tg ~ 1 n - 41, « 0 = 4-4 x 10 -8 X He 2 , 

Tg ~ 2 n - 19, £„ = 15 X He 2 . 

The total helium burning energy generation rate is 


L = 4ttRi 3 6 o (4.67) 

where 

J n = J*o x2 ^ 1_x ) n ^ (l+2x-1.8x 2 ) n dx, 

since 

T(r) - T c [1 + x - 3.8 x 2 + 1.8 x 3 ] . 


Then 


,t" ■ 201 6(jJ. 

L o 




The luminosity, with electron scattering opacity is igiven 


£ 


\v ) 
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by equation (3*16). The central temperature and density 
are from equations (3*6) and (3*3) 

.6 


•“ (i) ft) 1 • 


9.62 x 10 

5 


8 


The core structure, for T c ~ 1 x 10 °K, is 

8 „ 0*213 (g) 

„ -2.36 /MA 

u ° vv 


1.16 x 10 |i 


9.9 x 10' 


0.128 


- 1.62 


R 1 


8.27 x 10-V 0 - 787 /M ^ 0,872 


( 4 . 69 ) 


= 179 u 4 

L 0 


/ Ml\ 

\* o ) 


(SA 

I.V 


The density and temperature at the shell are determined 

by the conditions of hydrostatic equilibrium and energy 

conservation. The shell temperature is given by equation 

( 4 . 45 ) with the small second term neglected and equation 

( 4 . 69 ) for R, , _ - 

1 ^ xO .128 


T x = 6.97 x 10 


7 <^c°- 787 >d) 


. ( 4 . 70 ) 


The shell density is the solution of 

Luminosity - L core + L s hell> (4.71) 

where L core is the core luminosity, equation ( 4 . 69 ) > ^shell 
is the. shell energy generation rate, equation (4*50), and 
the luminosity is given by equations (4*51 ) » (3*24) or (3*25). 
Once is known^the luminosity is found from equation (4.71). 
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The effective temperature is given by equation ( 3 . 17 ) 

and the radius is given by equation ( 4 . 49 ) as in the case of 

hydrogen shell burning. 

When helium burning commences in the core of a star, 
the core expands, the central density . decreases, and the 
envelope contracts (Figure 16) . In massive stars, where 
the core was not degenerate, this adjustment is slight. In 
small mass stars, which developed a degenerate core during 
the helium core contraction, a helium flash occurs in which 

the core becomes nondegenerate and the central density is 

/ 

greatly reduced, the core expands and the envelope contracts. 
The resulting radii are much smaller than in the red giant 
stage, but still much larger than when on the main sequence. 

For small mass stars that have passed through the 
red giant stage, the luminosity during central helium 
burning is insensitive to the mass. The luminosity depends 
only on the core mass Usee equations (4 . 69 ) and ( 4 . 70 ) 3 , which 
is approximately the same at the onset of the helium burning 
in all such stars, since the smaller the mass of the star 
the larger the fraction of mass in the core. Thus small 
mass stars lie at. the onset of central helium burning in 
a strip of nearly constant luminosity but with varying 
effective temperature depending on the mass. 

The locus of points in the H-R diagram where initial 
central helium burning occurs is shown as the shaded 
regions in Figures 17 to 19. The relative contributions 
of hydrogen and helium burning to the luminosity are 
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found to be: Ljj » Ljj e for population X stars, while 

Ljj « Ljj e for population II stars. The time scales for 
evolution during the central helium .burning stage are 
given in Table 1. 

In the more advanced stages of evolution— helium 
burning, carbon burning, neon and oxygen burning — the 
core of the star continues becoming denser and hotter, 
a complicated shell structure develops, with some shells 
active and others inactive, and the radius continues to 
grow. A schematic picture of the stages of centrals 
nuclear burning and shell formation is given in Figure 

4 1 

5 of C. Hayashi, "Advanced Stages of Evolution," this 
conference, p. . How far a star progresses through 
these stages of nuclear burning depends, as we have shown, 
on its mass. 
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V. FINAL STAGES OF EVOLUTION 

After a star has exhausted all the nuclear fuels it 
is capable of burning, its only remaining sources of 
energy are its gravitational potential energy, which it 
can release by contracting, and its thermal energy, which 
it can release by cooling. Such a star will contract, 
increasing its central density and temperature. The core 
will, however, tend tol.be cooled off by energy losses from 
neutrino emission. The rate of emission of neutrinos 
increases with temperature, and since their mean free path 
is larger than the radius of the star they remove energy 
. from the star. If neutrino pair emission is intense, 
all stars in the stage of gravitational contraction after ' 
the exhaustion of nuclear fuel will develop degenerate 
cores. 

If the central density resulting from "the gravitational 
contraction is. low, only electrons, not nucleons, are 
degenerate and supply the pressure to support the star. 

There is a maximum density .possible for a stable star 
supported by degenerate electron pressure. At higher 
densities the electrons are forced onto the protons, 
creating neutrons. This process is a phase change and 
absorbs a great. deal of energy, causing instability. 

The gravitational collapse of massive stars produces 
cores with densities above the critical density. The core 
of such a star will be composed of free degenerate neutrons 
and other baryons. If the mass of the remnant from the 
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gravitational collapse is small enough, it can be sup- 
ported by the pressure of the degenerate neutrons and 
a stable neutron star will be formed. If the mass is 
too large, the gravitational force, augmented by the 

X 

relativistic effect that the pressure contributes to the 

j . 

effective mass, overwhelms the. nuclear forces, and the 

) 

star collapses indefinitely. What happens to such core 
remnants remains to be discovered. 

i 

\ 

Structure of White Dwarfs 

White dwarfs are stars whose support is provided by 
the pressure of degenerate electrons throughout most of 
. the mass of the star. In white dwarfs only electrons, 
not nucleons, are degenerate. We assume that the electrons 
are completely degenerate. This is, of course, not possible 
since in the surface' layers the density is very low and 
the electrons are nondegenerate. However, the surface 
layers are extremely thin. 

The equation of state of a degenerate gas is a compli- 
cated function 

P - P (p) 

approaching the limiting forms 

P - K x p 5/3 - 9.91 x 10 12 (p/n e ) 5/3 (5.1) 

at low density where the electrons are nonrelativistic 
(& « m e c), and 

P - K 2 p 4/3 * 1.23 x 10 1S (p/u e > 4/3 (5.2) 
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the mass, so the two forces will be in balance for only 
one mass, the limiting mass of a white dwarf star. For 
larger masses the gravitational force always exceeds the 
pressure force. 

The mass-radius relation for a white dwarf can be 
obtained from the virial theorem, (Salpeter, 1964) 


3 (y-l) U + fl =0} (2.1*) 

where Q is the gravitational potential energy, given by 
equation (1.19), 




GM 2 

R J 


and the internal energy U is the electron kinetic energy 

U = N K e 


Here N is the number of electrons and K is the kinetic 

e 

energy per electron. The mass of the star is 

M = N U e m p , 


where U is the molecular weight per electron 
© 



[x xi +| (y 1 + 2y 2 ) + |3 _1 , 


so — 2 for a fullv ionized gas if X = 0. Here m_ is 

e • p 

the proton mass. Thus from the virial theorem 



Gm p 2 u e 2 N 

3 (y-1)R 


( 5 . 4 ) 


The electron kinetic energy is related to its momentum 
by 

P e « m e c , 


K e = 


2m, 


K e * Pe c 


• P e >:> m e c 
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The average electron momentum p e is related to the average 
interelectron spacing r e by the uncertainty principle 

** e P e * h • 

Using the equality sign gives for the kinetic energy 

2 i r *-\2 


K, 


2 

£ 


Pe = i 
2m 


= 2 m e c 


K e = P e c 


m e c 


(?;) Pe <<C m e c l 

N (5.5) 

(r?) [?e >> m e4 


where r Q = ft/ m e c is the electron Compton wavelength. 

These two limiting equations can be combined in the inter- 
polation formula (Wheeler, 1964 ) 


K e " m e c ' 


by • 


(5.6) 


where s = r e / 


e/‘o * 


This formula is accurate to within 8/6. 


The radius of the star is expressed in terms of r by 

c 


R - N 1//3 r. 


(5.7) 


Equating the expressions (5*4) and (5.6) for the electron 
kinetic energy gives the relation 

,-2/3 

“P "* e " (5.8) 


1 + 2s 


f 0 3(Y-l)r 0 m e c^ 

“ — Z — F— N 

L Gm p u e -1 

- 3(y-1) f»o\ 2 '' 3 - 3 (y-1) I M cA 2 / 3 

u * '• N ' V 1 / ’ 


where 


■ m 


-3/2 


2.2 X 10 57 , 


\ 

$ 
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/ n x-3/2 

M o S («o) V = N o m p - 3.7xlo33 g 

* f 

i 

I 

- 1.85 % . ’ 

r 

For a nonrelativistic electron gas Y “ 5/3, and for an 
extreme relativistic electron gas y = 4 / 3 . The variation 


of y is given by Schatzman (1958), 

. , (*) 2/3 + 

Y-l 


2s 


3 7 ^_ f/3 i ' 

C 4 ) + S i 

3 (Y-l) varies between 1 and 2 and the above expression 1 
can be replaced by 

3<y-d - Vr 2 > 

1 + S 

with a maximum error of 27 % • Thus 


1 + s 


v 


Ho\ 


2/3 


(5.8') 


First note that the minimum value of the left-hand 
side of the above relation is 1, so there is a maximum 
mass for a white dwarf 


M 


max 


u - 2 M - 1 - 85 > 


(5.9) 


However, long before the density becomes infinite, inverse 
0 reactions will occur and the above analysis will cease 
to apply. The increasing density causes instability of 
the white dwarf before the singularity is reached. 

Second, the above relation can be written as a mass- 
radius relation 

-1/3 r -4/3 /M y 1 ^ /m \i/3 


R -■ R u 
o 


Cu 


e 


&) -&n* 


where 


1/3 8 

R Q ** N c r 0 = 5 x 10 cm . 

Thus the radius of a white dwarf is very small and it 
decreases as the mass increases. ' 

The mean density of a white dwarf is 



M/C-f- R 3 ) 


3MoUe r., -4/3 -3 


4TT R, 


1 


Gfe) 


- 1] 


7.06 x 10 u e Cu e 


-4/3 /M 


(fc) 


-2/3 


- 1] 


-3 


Since a white dwarf has a very thin nondegenerate surface 
layer, we may approximate it by a homogeneous model with 
a linear density distribution. Then the central density 
is ^ 

p c = 4P = 2.83 x 10 7 u e [H e ’ 4/3 - l] " 3 . (5.11) 

There is a maximum density possible for a stable white 
dwarf. As the density increases the electron Fermi energy 

increases. An electron with energy greater than the 

•/ ’■ ' 

p-decay -energy ,f or electron emission from a nucleus’ (Z-l, A) 

i ; i i | :1 

will produce inverse 0-reactions 1 

e - + (Z,A) - (Z-l,, A) + v . 


This process increases the value of u e in the interior, 

and thus the maximum stable mass is reduced. The predominant 

nuclei under white dwarf conditions are elements in the 
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!'• 

i ; 

I ' V 

i 

range neon to iron, for which inverse 3-decay will 
occur at densities about 10^ g/cm^ . Thus the critical 
density for a white dwarf is about 10^ g/cm^ . The relation 
between central density and mass for a white dwarf is 
shown in Figure 22 , from Wheeler (1964) . The stable con- 
figurations shown at higher densities are the neutron 
stars . 

The degenerate interior of a white dwarf is practically 

isothermal because heat conduction by degenerate electrons 

1 

is very efficient. This isothermal interior is blanketed 
by a nondegenerate surface layer, which is very thin and 
contains only a minute fraction of the mass of the star. 

The small extent of the surface layer is easily seen by 

t 

considering the scale height 

P RT 

Pg dg- 

The temperature at the transition layer is df the order 
of a million degrees^ but g = GM/R^ is extremely large \ 
because R is very small. Assuming M « Mq, R « R q and 
T » 10 6 , then g » 5 x 10 8 and t » 10 6 cm = 10 km. The 
density in the surface layer is less than about 10^ g/cm^ 
since it is nondegenerate. Again assuming T 10 ^, the 
mass of the surface layer will be 

M s = 4TTR 2 p AR ss tt 10 27 10“ 6 . 

Therefore the equations for the surface layers may be 
integrated explicitly since g, M, and L are practically 


constant 



I 


Figure 22 . Schematic mass density relation for white dwarfs 

and denser configurations from Wheeler (1964), calculated 
for cold catalyzed matter. 


\ 
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The surface layer is in hydrostatic equilibrium and 
energy transport is by radiation.. We will assume Kramer's 


opacity. 


-3.5 


H = K Q T 
o 

K q = 4-34 X 10 


25 Z (l+X) 

(t/g) ' 


where (t/g) is a quantum mechanical correction factor w 10 
in this case. Then the equations for the structure of the 
envelope are (Schwarzschild, 1958 and Chandrasekhar, 1939) 


72 P 


16d 4 ff r^ * 

64 rraGM ^3 
3 xL 


= 64 ttctGM k p -l t 7.5 # 

3k q L uH 

Thus the pressure and density are related to the temperature 


2 64 rroGMk " \ 2 4.25 

.5 3K 0 IUH j 


( 5 . 12 ) 


" (K 


2 64 ttctGMUH ) 2 ,, t,3 .25 

13T5 3'VoLk "7 T 


The radial dependence of T can be found from the equation 
of hydrostatic equilibrium 


dP = 
dr 


p UH 1 

f IT GM 72 


and from equation ( 5 * 12 ) 


4.25 


Thus 
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(5,-13 ) 


These equations for T, P, and p can be used throughout 

j 

the nondegenerate surface layer. 

The properties of the transition layer between the 
degenerate interior and the nondegenerate surface layer 
can be found as a function of the luminosity of the white 
dwarf (Schwarzschild, 1958). The isothermal nature of 
the interior gives a relation between interior temperature 
and the luminosity, which is constant through the surface 
layers, as follows J 


L = - 2 


875 


64 ttqGM uH 
3 k 


r>6 « 5 


(5.14) 


Apply this to the transition layer. The boundary condition 
is the equality of the electron pressures in the two regions 



so the boundary condition is 



2.4 x 10 ^ u e 



(5.15) 


Then the luminosity and internal temperature, T_ = 

c 

related by 


L 


2 64 TTCfGM /H \ 4 u T 3,5 

8 - 5 "V k 0 k7 UJ 7? c 

5.7 = 10 25 (^)-K T ,^ 3 - 5 

U / Ue c (6) 


T tr , are 


( 5 . 16 ) 
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The internal temperature, transition density, and extent 
of surface layer as a function of luminosity is shown in 
the following table for a 1 solar mass star with composition 
X - 0, Y = 0.9, Z - 0.1. 


1 


log Ptr 

R - r tr 
R 

cs 

* 

o 

H 

■ 

17 

i 

3.5 

0.011 

10" 3 

• 9 

j 311 

0.006 

1 

O 

H 

4 

! 2.6 

t 

0.003 


This table is taken from M. Schwa rzschi Id: Structure 

and Evolution of Stars . - p. 238. 

‘ t ■ 

* * * 

The source of energy for white dwarfs is the thermal i ji 
energy of the nondegenerate nuclei. The energy source 
cannot be nuclear reactions. At the high densities found 
in white dwarf interiors the Coulomb barriers of nuclei 
are reduced. At densities greater than about 5 ^ 10 g/ cm 
hydrogen reactions occur and at densities greater than 
about 5 x 10^ g/cm^ helium reactions occur. However, 
during a star 1 s evolution before becoming a white dwarf, 
all the hydrogen in its core will have been exhausted, 
while white dwarfs with central densities high enough for 

heliUm reactions are massive enough to have exhausted the 

\ ' 

helium in their cores. In the surface layer, where hydrogen 
may be abundant, nuclear reactions would cause instability 
because of their temperature sensitivity. During a con- 
traction, the rate of energy generation would increase 
above and during an expansion, would decrease below Its 
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equilibrium value, thus feeding energy into the pulsations. 
The energy source cannot be gravitational, because a star's 
radius is fixed by the mass-radius relation after it has 
become .almost completely degenerate and no further con- 
traction is possible. The energy source cannot be the 
thermal energy of the electrons because they are degenerate 
and most are already in their lowest possible energy state. 

The evolution of a white dwarf is a continual slow 
cooling at constant radius; its luminosity and effective 
temperature decrease in time. Evolutionary paths in 


the H-R diagram are shown for several masses in Figure 23. 

I 

I^e luminosity of a white dwarf is the rate of change 
of the thermal energy of the nondegenerate nuclei. 



d 

dt 




( 5 . 17 ) 


where is the molecular weight of the nuclei, 

-1 _ j 

— X + ^Y. This equation can be integrated to obtain 
the cooling time of a white dwarf (Schwarzschild, 1958) . 
Using the expression for the luminosity, (5.16)> 

L - K(n,M) T 3 * 5 


gives 


dT 

dt 


C T n , 


where n — 3*5 and C — - ^ * Integration gives the 

cooling time from "infinite" temperature, setting the 

4 

integration constant equal to zero, which is the time 
scale of evolution of white dwarfs 



i 

! 


i 

i 

* ' 1 

i 

Figure 23. Evolutionary tracks of white dwarfs in the H-R 

diagram. Solid curves are from analytic expression (5.10), 
the dashed curve is from Schwarz schi Id, 1953 . ■ 
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